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Abstract 

In shallow-water waveguides a propagating field can be decomposed over three kinds of modes: 
the propagating modes, the radiating modes and the evanescent modes. In this paper we consider 
the propagation of a wave in a randomly perturbed waveguide and we analyze the coupling 
between these three kinds of modes using an asymptotic analysis based on a separation of scales 
technique. Then, we derive the asymptotic form of the distribution of the mode amplitudes 
and the coupled power equation for propagating modes. From this equation, we show that the 
total energy carried by the propagating modes decreases exponentially with the size of the random 
section and we give an expression of the decay rate. Moreover, we show that the mean propagating 
mode powers converge to the solution of a diffusion equation in the high-frequency regime. 
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Introduction. 

Acoustic wave propagation in shallow-water waveguides has been studied for a long time because of 
its numerous domains of applications. One of the most important applications is submarine detection 
with active or passive sonars, but it can also be used in underwater communication, mines or ar- 
chaeological artifacts detection, and to study the ocean's structure or ocean biology. Shallow-waters 
are complicated media because they have indices of refraction with spatial and time dependences. 
However, the sound speed in water, which is about 1500 m/s, is sufficiently large with respect to the 
motions of water masses that we can consider this medium as being time independent. Moreover, 
the presence of spatial inhomogencities in the water produces a mode coupling which can induce 
significant effects over large propagation distances. 

In shallow-water waveguides the transverse section can be represented as a semi-infinite interval 
(see Figure [lj and then a wave field can be decomposed over three kinds of modes: the propagating 
modes which propagate over long distances, the evanescent modes which decrease exponentially with 
the propagation distance, and the radiating modes representing modes which penetrate under the 
bottom of the water. The main purpose of this paper is to analyze how the propagating mode powers 
are affected by the radiating and evanescent modes. This analysis is carried out using an asymptotic 
analysis based on a separation of scale technique, where the wavelength and the correlation lengths 
of the inhomogencities, which are of the same order, are small compared to the propagation distance, 
and the fluctuations of the medium are small compared to the wavelength. In the terminology of |7J 
this is the so-called weakly heterogeneous regime. 

Wave propagation in random waveguides with a bounded cross-section and Dirichlet boundary 
conditions (see Figure [T]) has been studied in [7j Chapter 20] or [9] for instance. In this case we have 
only two kinds of modes, the propagating and the evanescent modes. In such a model an asymptotic 
analysis of the mode powers show total energy conservation and an equipartition of the energy carried 
by the propagating modes. In [5] coupled power equations are derived under the assumption that 
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Figure 1: Illustration of two kinds of waveguides. In (a) we represent a shallow- water waveguide 
model with an unbounded cross-section. In (6) we represent a waveguide with a bounded cross-section. 



evanescent modes are negligible. In |5] the role of evanescent modes is studied in the absence of 
radiating modes. In this paper we take into account the influence of the radiating and the evanescent 
modes on the coupled power equations. In this case we show a mode-dependent and frequency- 
dependent attenuation on the propagating modes in Theorem |5.1| that is, the total energy carried by 
the propagating modes decreases exponentially with the size of the random section and we give an 
expression of the decay rate. Moreover, in the high-frequency regime, we show in Theorems |5.2| and 



5.4 that the propagating mode powers converge to the solution of a diffusion equation. All the results 



of this paper are also valid for electromagnetic wave propagation in dielectric waveguides and optical 
fibers [TZl HH1 EH EH [29] . 

The organization of this paper is as follows. In Section [l] we present the waveguide model, and 
in Section [2] we present the mode decomposition associated to that model and studied in detail in 
|28j . In Section [3] we study the mode coupling when the three kinds of modes are taken into account. 
In the same spirit as in [7J Chapter 20], we derive the coupled mode equations, we study the energy 
flux for the propagating and the radiating modes, and the influence of the evanescent modes on the 
two other kinds of modes. In Section [4] under the forward scattering approximation, we study the 
asymptotic form of the joint distribution of the propagating and radiating mode amplitudes. We apply 
this result in Section [5] to derive the coupled power equations for the propagating modes, which was 
already obtained in [2] or [TH] for instance. In this section, we study the influence of the radiating 
and evanescent modes on the mean propagating mode powers. We show that the total energy carried 
by the propagating modes decreases exponentially with the size of the random section and we give 
an expression of the decay rate. In other words, the radiating modes induce a mode-dependent 
attenuation on the propagating modes, that is why these modes are sometimes called dissipative 
modes. Moreover, under the assumption that nearest-neighbor coupling is the main power transfer 
mechanism, we show, in the high-frequency regime or in the limit of large number of propagating 
modes, that the mean propagating mode powers converge to the solution of a diffusion equation. We 
can refer to [21 [TH] for further references and discussions about diffusion models. In that regime, we 
can also observe the exponential decay behavior caused by the radiative loss. 



1 Waveguide Model 

We consider a two-dimensional linear acoustic wave model. The conservation equations of mass and 
linear momentum are given by 



, du „ „ 

p(x,z)—+X7p = F 

1 dp „ 

1 V.u = 0, 



(1) 



K(x, z) dt 



where p is the acoustic pressure, u is the acoustic velocity, p is the density of the medium, K is 
the bulk modulus, and the source is modeled by the forcing term F(t,x,z). The third coordinate z 
represents the propagation axis along the waveguide. The transverse section of the waveguide is the 
semi-infinite interval [0,+oo), and x S [0,+oo) represents the transverse coordinate. Let d > 0, we 



C. GOMEZ 



SHALLOW-WATER PROPAGATION 3 




Figure 2: Illustration of the shallow-water waveguide model. 



assume that the medium parameters are given by 



!±(n 2 {x) + *JeV{x,z)) if x€[0,d], z€[0,L/e] 
( x S [0,+oo), z £ (-oo,0) U (L/e,+oo) 
^n 2 (x) if I or 

| i£ (d, +oo), z £ (— oo,+oo). 

p(x, z) — p if x £ [0, +oo), z£l. 

In this paper we consider the Pekeris waveguide model. This kind of model has been studied for half 
a century [23 and in this model the index of refraction n(x) is given by 



n{x) 



n\ > 1 if x £ [0, d) 
1 if ie [d, +oo). 



This profile can model an ocean with a constant sound speed. Such conditions can be found 
during the winter in Earth's mid latitudes and in water shallower than about 30 meters. The Pekeris 
profile leads us to simplified algebra but it underestimates the complexity of the medium. However, 
the analysis that we present in this paper can be extended to more general profiles n(x) with general 
boundary conditions. In the Pekeris model that we consider n\ represents the index of refraction of 
the ocean section [0, d], where d is the depth of the ocean, and we consider that the index of refraction 
of the bottom of the ocean is equal to 1. This model can also be used to study the propagation 
of electromagnetic waves in a dielectric slab and an optical fiber with randomly perturbed index of 
refraction [HJ HH • 

We consider a source that emits a signal in the z-direction, which is localized in the plane z = L$. 



F(t,x,z) = V(t,x)6(z- L s )e z 



(2) 



fy(t, x) represents the profile of the source and e z is the unit vector pointing in the ^-direction. < 
is the location of the source on the propagating axis. 

The random process (V(x, z),x £ [0,d],z > 0) that we consider, and which represents the spatial 
inhomogeneities is presented in Section [6. 1| However, one can remark that the process V is unbounded. 
This fact implies that the bulk modulus can take negative values. In order to avoid this situation, we 
can work on the event 



(v(a, z) £ [0, d] x [0, L/e], m + VtV(x, z) > o). 



In fact, the property (55 1 implies 



lim F(3(x, z) £ [0, d] x [0, L/e] : m + y/eV{x, z) < o) 



< lim 



imP(Ve 



sup sup 

ze[0,L] a;e[0,d] 



V 



K) 



> n- 
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2 Wave Propagation in a Homogeneous Waveguide 

In this section, we assume that the medium parameters are given by 

p(x, z) — p and K(x, z) = , V(x, z) € [0, +oo) x M. 

n (x) 

From the conservation equations ([TJ , we can derive the wave equation for the pressure field, 



where c(x) = c/n{x) with c = y y , an d A = c* 2 . + d 2 

In underwater acoustics the density of air is very small compared to the density of water, then it 
is natural to use a pressure-release condition. The pressure is very weak outside the waveguide, and 
by continuity, the pressure is zero at the free surface x = 0. This consideration leads us to consider 
the Dirichlet boundary conditions 

p{t,0,z) = V(i, z) e [0, +oo) x E. 

Throughout this manuscript, we consider linear models of propagation. Therefore, the pressure 
p(t, x, z) can be expressed as the superposition of monochromatic waves by taking its Fourier trans- 
form. Here, the Fourier transform and the inverse Fourier transform, with respect to time, are defined 
by 

/M = J f(t)e^dt, f(t) = i- J fi^e-^dw. 

In the half-space z > L$ (resp., z < L$), taking the Fourier transform in (ujj, we get that p(u>, x, z) 
satisfies the time-harmonic wave equation without source term 

d 2 z p{oj, x, z) + d 2 x p{uj, x, z) + k 2 (uj)n 2 (x)p(uj, x, z) = 0, (4) 

where fc(w) = — is the wavenumber, and with Dirichlet boundary conditions p(u>, 0, z) — Vz. The 
source term implies the following jump conditions for the pressure field across the plane z = L$ 

p(u,x,L^) -p(uj,x,L s ) = $(u,x), 
d z p(uJ,x,L^) - d z p(u,x,Lg) = 0. 

2.1 Spectral Decomposition in Unperturbed Waveguides 

This section is devoted to the presentation of the spectral decomposition of the Pekeris operator 
d 2 + k 2 (ui)n 2 (x) . The spectral analysis of this operator is carried out in [28]. Throughout this paper 
we are interested in solutions of Q such that 

p(cj, ., .)l(£ s ,+oo)(z) 6 C° ((£<?, +oo),fl 1 (0, +c^) n H 2 (0, +oo)) n C 2 ((L 5 , +oo), Hj , 
p(u, ., .)l ( _ co , Ls) (z) G C°((-oo, L s ) t H^0, +oo) n H 2 (0, +oo)) n C 2 ((-c», L s ), h) , 

where H = L 2 (0, +oo). H is equipped with the inner product defined by 

V(fti, h 2 ) <E H x H, (hi,h 2 ) H = / h 1 (x)h 2 {x)dx. 

Jo 



Consequently, in the half-space z > L$ (resp., z < Ls), we can consider Q as the operational 
differential equation 

d 2 

^p{uj,.,z) + R(uj)(p{u>,.,z))=0 (6) 
in H, where R(w) is an unbounded operator on H with domain 

V(R{u>)) = H%(0, +oo)nH 2 (0,+oo), 
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and defined by 



R(Lu)(y) = ^V + k 2 {uo)n 2 {x)y Vy e V(R(u)). 



According to [25], R(ui) is a self-adjoint operator on the Hilbert space H, and its spectrum is given 
by 

Sp(R(u)) = (-«3,fc 2 M] U{^ (w) (a;),... ,/?!»}. (7) 
More precisely Wj € {l, ... , A^w)}, the modal wavenumber /3j(w)is positive and 

fc 2 M < < • • • < PKlo) < n\k 2 (u). 

Moreover, there exists a resolution of the identity II U of R(oj) such that Vy £ -ff, Vr £ K, 



AT(w) 



IL(r,+oo)(y)(x) = (^iC^O^jK^r.+cc) (AH 2 ) 
J'=l 

{yA~t{u, .)) i? ^ 7 (w,x)d7l(_ 00ife 2 (w ))(r), 



and Vy € V(R(u)), Vr € R, 



JV(w) 



n w (r,+oo)(i?M(y))0) = 51 /M w ) 2 (j/>0j( w >O)#<M w > x ) 1 (r,+oo) (ftM ; 



Let us describe these decompositions. 



Discrete part of the decomposition Vj ' E {!,..., N(cu)}, the jth eigenvector is given by 



Aj(u>)sm.(<jj((jj)x/d) if < x < d 

^■(w)sin(CT J (u;))e"^( a) )^ if d < x, 



where 
and 



2/d 



A . sin 2 (o"j (cj)) sin(2crj (cj)) 

According to [28 , <7i(a;), . . . ,cr/v( w )(w) are the solutions on (0, nik(uS)dQ) of the equation 

y 



tan(y) = - 



V(niW-y2' 



(8) 



(9) 



such that < cri(tj) < • • • < ct n ^(lu) < nik(Lu)dd, and with 9 = \J\ — l/n\. This last equation ad- 
mits exactly one solution over each interval of the form (jr/2+(j — 1)tt, 7r/2+j7r) for jg {1, . . . , A(w)}, 
where 

7T 

and [•] stands for the integer part. From ([9]), we get the following results about the localization of the 
solutions which is used to show the main result of Section 15.21 

Lemma 2.1 Let a > 1/3, we have as N(ui) — > +00 



sup 

jE{l,...,N(u)-[N(u,)<*]-l} 



sup 

j£{l,...,JV(w)-[JV(a.)«]-2} 



\a j+1 {u)-a j {u)-Tr\ = 0(N{u)?-i°). 
<r j+2 (u) - 2<r i+ i( W ) + <r»)| - O (iV^) 1 - 3 ") 
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Continuous part of the decomposition For 7 € (— 00, k 2 (u>)), we have [25] 

(/)y(u,x) = 

{Ay(uj)sin(r](u!)x/d) if < x < d 

A 7 H(sin(r ? H)cos(eH^) + f^cos(^H)sin(eH^)) if d < x, 

where 

r?(w) = d^Jn 2 k 2 (u) - 7, £(cj) = dy/k 2 (uj) - 7, 

and 



Y ^(e 2 Hsin 2 (7 ? H)+7 7 2 Hcos 2 (7 ? H))' (10) 
It is easy to check that the function 7 1— > A 7 (w) is continuous on ( — 00, k 2 (oj)) and 



AH ~ 7=rn7i - ( n ) 



We can remark that <f>y(u>, .) does not belong to H . Then, (y, <f>^(u>, -)) H is not defined in the classical 
way. In fact, 

(y,4> 1 {uj, .)) „ = lim / y(x)(f> 1 (u),x)dx on L 2 ( — 00, k 2 (cu)) . 
Moreover, we have Vy £ H 

N X ] 2 /" fc2(Lj) 2 

l|y||ff= 2^ Ky><M w >-)} H | + / \{vA-y{^,-)) H \ ^7- 

Then, 

6 W : i? — ► ft" 

y — » (((^^^^hIh ^.(^^(^-llfll^f-oo^M)) 

is an isometry, from H onto ft" = xi 2 (-oo, k 2 (u)). 

2.2 Modal Decomposition 

In this section we apply the spectral decomposition introduced in Section [2d] on a solution p(u>, x, z) 
of the equation (|6]). Consequently, we get the modal decomposition for p(uj,x,z) in the half-space 

z > L s , 

p\u,x,z) = ^2p j (u,z)(f> j (uj,x)+ / p 7 (u,z)(f>~ / (cj,x)d'y, 

3 = 1 J -°° 

where p(uj, z) = 6^(^(0;, ., z)). For j £ {l, . . . , N(cu)}, S^oII^dj}) represents the projection onto the 
jth propagating mode, and Pj(u>, z) is the amplitude of the jth propagating mode. Q u o II w (0, k 2 (uj)) 
represents the projection onto the radiating modes, and p 7 (w, z) is the amplitude of the 7th radiating 
mode for almost every 7 £ (0,k 2 (uj)). Finally, O w o U u (— 00, 0) represents the projection onto the 
evanescent modes and p 7 (cj,z) is the amplitude of the 7th evanescent mode for almost every 7 £ 
(-oo,0). 

Consequently, p(u, z) satisfies 

d 2 

^Pj(w, z) + Pj (Lv)pj(uj 7 z) = 0, 
d 2 

j-^P 7 (u, z) + 7 p 7 (w, z) = 
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in 7i w and the pressure field can be written as an expansion over the complete set of modes 



^p fe iftMz J ^,x) + 



£7,0 M 



(j>j(uj, x)dj 



L (is,+oo) 



~[ V&M 7o 7 i/4 



(12) 



L(-oo,i 



; )( Z )' 



under the assumption that the coefficients (c 7i o(w)e v ^' 7 ' is /|7| 1 ^ 4 ) 7 and (d^^(uj)e"y^ Ls /l^l 1 ^ 4 ^ 
belong to L 2 (—oo, 0). 

In the previous decomposition, aj o(ui) (resp., bjfi(u)) is the amplitude of the jth right-going (resp., 
left-going) mode propagating in the right half-space z > L$ (resp., left half-space z < Lg), a 7] o(w) 
(resp., 6 7i o(w)) is the amplitude of the 7th right-going (resp., left-going) mode radiating in the right 
half-space z > L$ (resp., left half-space z < L$), and c 7) o(w) (resp., djfl(uj)) is the amplitude of the 
7th right-going (resp., left-going) evanescent mode in the right half-space z > L5 (resp., left half-space 
z < L s ). 

We assume that the profile ^(t,x) of the source term Q is given, in the frequency domain, by 



■9(u,x) = f(oj) 



3=1 



(-S,-«)U(«J,fc 2 (o;)) 



<pJu,xo)<f>-y(uj,x)d~/ 



(13) 



where xq £ (0, cf). The bound S in the spectral decomposition of the source profile was introduced to 
have ^(u>, .) E H, and £ was introduced for technical reasons. Note that S can be arbitrarily large 
and £ can be arbitrarily small. Therefore, the spatial profile in ( 13 1 is an approximation of a Dirac 
distribution at xq, which models a point source at xq. 
Applying 6 W on ^ and using ( [T2| , we get 



a^H - -&;,oM = ^^7(c)^(c,x )e-^(^ Ls Vj e {l, . . 

a (u) = -6 o(lj) = I 2 T-/( w )^7( w ^o)e" 4 ^ is for almost every 7 G (£, k 2 {u)) 
1 for almost every 76(0, £), 

c 7 ,oM = ^/>)^( Wl3 ;o)e^ is , d 7 , H = 2^/( w )^ 7 ( w , xo)e~^ is 
for almost every 7 S (— S, — £), and 

c 7)0 (w) = (iy.oM = 
for almost every 7 £ (—00, —5) U (— £, 0). 



3 Mode Coupling in Random Waveguides 

In this section we study the expansion of p(lo, x, z) when a random section [0, L/e] is inserted between 
two homogeneous waveguides (see Figure [2]) . In this section the medium parameters are given by 

!j ? (n 2 (x) + ^eV(x,z)) if xe[Q,d], z€[0,L/e] 
( x e [0,+oo), z e (-oo,0) U (L/e, +00) 
in 2 (a;) if I or 

I i£ (d, +00), z 6 (—00, +00). 

p(x, z) — p if x £ [0, +00), zel. 
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In the perturbed section, the pressure field can be decomposed using the resolution of the identity 
II W of the unperturbed waveguide: 

p{w,x,z) = y ^2,p j (u],z)<j)j{u,x)+ I p 7 (w,z)0 7 (w,a;)d7, 

3 = 1 J -°° 

where p(cu, z) = w (p(w, ., z)). In what follows, we shall consider solutions of the form 

p(uj,x,z)= y2pj(u,z)<j)j(u,x)+ / p 7 (uj,z)(j) 7 (uj 1 x)d'y. 
j = l J(-oo,-0u(C,fc 2 (^)) 

This assumption leads us to simplified algebra in the proof of Theorem |4.1| In such a decomposition, 
the radiating and the evanescent part are separated by the small band (—£,£) with (< 1. The goal 
is to isolate the transition mode between the radiating and the evanescent part of the spectrum 
Sp{R(ij)) given by d7L Moreover, we assume that e <C £ and therefore we have two distinct scales. 
Let us remark that in this paper, we consider in a first time the asymptotic e goes to and in a second 
time the asymptotic £ goes to 0. 

3.1 Coupled Mode Equations 

In this section we give the coupled mode equations, which describe the coupling mechanism between 
the amplitudes of the three kinds of modes. 

In the random section [0,L/e] the pressure field p(oj,z) satisfies the following coupled equations 
in W: 

d 2 N[UJ) 
j^Pj(u, z) + f3j((j)pj(u, z) + yfek 2 (uj) ^2 Cji{ z )Pl{u, z) 

i=i 



+ ^~ek 2 {u) / CV ,{ z )p ll {u,z)di = 0, 

J(-00 -?)U(€:fc 2 (u,)) 



2 JVM 



(14) 



d 

z) + 7 p 7 (w, z) + y/ek 2 (uj) ^ C^(z)pi(uj, z) 



dz 

1=1 



+ V~ek 2 (u;) / C",(z)py(w,z)dy = 0, 



where 



^l( z ) = (<i > 3(< J ,.),<t>i(w,-)V(.,z)) H = / <f>j(u,x)4>i(w,x)V(x,z)dx, 

Jo 

,d 

Cf{z) = C^j(z) = ((j>j{u,.),<j)^{u,.)V{.,z)) = / <j) j (uj,x)<j) 1 (uj,x)V(x,z)dx, (15) 

Jo 

f d 



We recall that p(u, .,.) £ C°((0, +oo), i^(0, +oo) n H 2 (0, +oo)) n C 2 ((0, +oo), H), then 

-€ 

7 2 |p 7 (w,z)| 2 ,i7 < +oo. (16) 

In the previous coupled equation the coefficients C u (z) represent the coupling between the three kinds 
of modes, which are the propagating, radiating and evanescent modes. 

Next, we introduce the amplitudes of the generalized right- and left-going modes a{uj, z) and 
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b{uj 1 z), which are given by 



dz 



%{u,, z) = iy/fwjfau, z)e^ z z)e- i ^ z ) , 



P 7 (w,z) = - T j i (cL 1 { tl ),z)e i ^" +b 7 (uj,z)e 
-p 7 (w,z) = ^7 1 / 4 (a 7 (^z)e ^ ^ -b^(u,z)e- i ^ 
Vj € {l, . . . , N(lo)} and almost every 7 e (£, k 2 (uj)). Let 



From p4| , we obtain the coupled mode equation in Ti^ x x L (— 00, — £) for the amplitudes 
(a(w, z), 6(^,2:), p(w,2:)): 



(=1 v/W 

i*» ^> ( Sy(WjZ)e <yy- ft ), + 6 7 ,( w ,, )e -^ + ft)^ H (17) 



> sfWi 



' e -i0jZ 



dz 



ifc 2 (a;) /- fc2(w) GSL,(«) 



2 / \ /^w 



(18) 



ik 2 (u;) f k2 ^ Cp(z) 



^fOZL (a Y (cu,z)e^+^> + b 1 ,{u J ,z)e-^-^) z ) dr/ 



(19) 



? fc 2 M r £ cv { Z ) 



p~f'(ui, z)d'y'e 



dz 2 i~i VV^Pi v ' 

- -fi^ f' M (a,( W ,*)e<(^). + 6 ¥ (»^)e- < (^-^>) */ < 20 > 

d 2 

— p 7 (uj, z) + 7 p 7 (w, z) + y/eg y (uj, z ) = °j ( 21 ) 
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where 



"K^j qui r z \ 

9l {w,z) = k 2 (cu) J2 {a l {uo,z)e^ + b l {o J ,z)e-^ z ) 

1 — i Vr/ 



<- fc V) CJ y (z) 
7 

+ fc 2 H/ C^,(z)p Y (u;,z)d 1 / . 



(22) 



DC 



Let us note that in absence of random perturbations, the amplitudes a(uj, z) and b(ui, z) are constant. 

We assume that a pulse is emitted at the source plane L$ < and propagates toward the randomly 
perturbed slab [0, L/e]. Using the previous section, the form of this incident field at z — is given by 

J=t V&H A 7 i/4 J-s l7l 1/4 

Consequently, by the continuity of the pressure field across the interfaces z = and z = L/e, the 
coupled mode system is complemented with the boundary conditions 

a(uj, 0) = ao(uj) and b ( u, — ] =0 



in Ti^. For j G {l, . . . , 7V(w)}, aj o(^) represents the initial amplitude of the jth propagating mode, 
and for 7 £ (£, k 2 (cu)), a 7 , (w) represents the initial amplitude of the 7th radiating mode at z = 0. 
Moreover, for 7 £ (— S 1 , — £), c 7i o(w) represents the initial amplitude of the 7th evanescent mode at 
z — 0. The second condition implies that no wave comes from the right homogeneous waveguide. 

3.2 Energy Flux for the Propagating and Radiating Modes 

In this section we study the energy flux for the propagating and radiating modes, and the influence 
of the evanescent modes on this flux. 

We begin this section by introducing the radiation condition for the evanescent modes 

lim |[n w (-oo,-£)(p(w,.,z))|| =0. 

This condition means, in the homogeneous right half-space, that the energy carried by the evanescent 
modes decay as the propagation distance becomes large. From the radiation condition and (21 1, we 
get for almost every 76 (— 00, — £) 

pJiu,z) = ^= f ' gJuj,u)e^ u ^ z Uu+ -^L f ' gJu,u)e^ z - u) du 



Vz £ [0,+oo). According to (12 1, the relation (24 1 can be viewed as a perturbation of the form of the 
evanescent mode without a random perturbation. Using the same arguments as in [7J Chapter 20], 
we get Vz £ [0, L/e], 

^;(l|a(^ z )ll^ - IMw.^Hw-) = -\felm ( j s 7 (w, z)p 7 (uj, z)dry J , 



and 



\\a(u;,z)\\ 2 -||6(c,z)|| 2 HlaoMH 2 -||6oM||^ / ? ^¥^7 



2 /wFfl _ (25) 

y/l I cj) 1 (u,x )eVW Ls [ Im(g 7 (uj,u))e-^ u dudj, 
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where 



G 7 (w, z) = 



z r L/e 



Jz 



Im(g 7 (u}, tt)g 7 (cj, «))e"v M(« v ^dvdu 



Consequently, for z = L/e, we get 



+||S( W ,0)||2 =||aoM|& 



5 " 

/-£ , /-i/e , 

(/> 7 (w,a; )evM L s / Im(g y (uj, u))e _ v ^ u c2m d 7 . 
-S ^0 



The second term on the right side of the previous relation has the factor 4> 1 {lo, Xo)e v ' 7 ' is which is the 
form of the evanescent mode at z = without a random perturbation. Therefore, if is is far away 
from and whatever the source (evanescent modes decay exponentially from L$ to 0) or if there is 
no excitation of modes 7 G (—00, — £) by the source (that is when S = £), we can get the conservation 
of the global energy flux for the propagating and radiating modes: 



\\a(o;,L/e)f H . + 116(^,0)11^ = \\M")\ 



However, from (251 and even if there is no evanescent modes in (23 1, the local energy flux is not 



conserved. The energy related to the evanescent modes is given by the last two terms on the right 



side in (p5|. Let us estimate these two quantities. First, 
• G~f(u>, z) 



sup 

*e[0,L/e] 



-c?7 



<K(£,d) sup sup v(x,-\ 

z€[0,L] x€[0 : d} ^ e ' 

x sup \\a(cj,z)\\^ + \\b(u,z)\\^ 

ze[0,L/e] ( 5 



IIF(^)I' 2 



Second, 



-£ 



7 (w,x o )eVl^l i s / Im( y g 1 ( U j,u))e-^^ u dud 1 



Vlx, 



sup 

zE[0,L/e] 

< K(£, d) sup sup 
ze[o,L] xe[o,d] 

x sup \\a(uj,z)\\ n ^ + \\b(u),z)\\nv + ||p(w, z)|| L i(_ 00 ._ 4) . 

ze[0,L/e] 

In the two previous inequalities if (£, d) represents a constant which can change between the different 
relations. However, it is difficult to get good a priori estimates about 



sup ||a(w,z)||L, + \\b(u),z)\\ n „ + \\p{u),z)\\ L i ( _, y 

ze[0,L/e] 5 5 



(26) 



For this reason, let us introduce the stopping "time" 



L e = m£[L>0, sup \\a(uj,z)\\^ + \\b(u,z)\\^ + WK^^Uhi-o^-c) 



> 



1 



ze[0,L/e] 

The role of this stopping "time" is to limit the size of the random section to ensure that the quantity 



(p6| is not too large. Consequently, the energy carried by the evanescent modes over the section 

2 



[0,L/e] for L < L c , is at most of order ©(e 1 / 4 sup z£ [ L / e j su Pa;e[o,d] 1^ ( x ' z )\ )> an(1 according to (551 
the local energy flux for the propagating and the radiating modes is conserved in the asymptotic 
e — > 0. More precisely, we can show that Vrj > 0, 



lim P sup 

\zG[0,L/e] 



\\a(w, z)\\* - \\b(u, z)f m - IISoMH^ + ||6oHfe > r), L < U = 0. (27) 



In Section |4j we shall see, under the forward scattering approximation, that the condition L < L e is 
readily fulfilled in the limit e — > 0, that is we have lim e _^o V(L e < L) = 0. 
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3.3 Influence of the Evanescent Modes on the Propagating and Radiating 
Modes 

We analyze, in this section, the influence of the evanescent modes on the coupling mechanism between 
the propagating and the radiating modes. 

First of all, we recall that 0^ o oo, — £) (p(u>, ., z)j represents the evanescent part of the 

pressure field p[u>,.,z), where 0^, and are defined in Section 2.1 In this section we consider 
F = L l {— oo, — £) equipped with the norm 



\y\\F 



\Vi\dl, 



which is a Banach space. Substituting (22 1 into (|24j), we get 

{Id - Ve$") (e w o n w (-oo, -0 (p(io, .,.))) = y/ip(w, .) + po(u, •)• 
This equation holds in the Banach space (C([0, +oo), F) , ||.||oo.f), where 

\\y\\oo <F = sup \\y(z)\\ F Vy G C([0, +oo), F) . 

z>Q 

In (281, $ w is a linear bounded operator, from (C([0, +oo), F) , ||.||oo.f) to itself, defined by 



(28) 



jfc 2 (w) 



zAL/e 



fc 2 H 



L/e 



2 a/ItI JzAL/eJ-ac 

Vz € [0,+oo), and for almost every 7 G (— 00,— £) 

fc2(w) t-zAL/e n-^j) C"j,(u) 



C^,(u)y y (u)dYevM(*-«Odu 



p 7 (w,z) 



Jo L r~7 VP; 



1 21 \ rL/e r N (") ^ , s 

" M ' [2^#(ai( W ,«)e**« + 6i( W ,«)e-^) 



2 V |-y| «A*a£/ 



( ) G 7T7?(«7'( w ^) eW " + 6 7 '(w,u)e- l ^ 7u )]d 7 'e^ (z - u) dii 



Vz G [0,+oo). Finally, for almost every 7 G (—00, — £) and Vz G [0,+oo), 

m.otw,*) = 7 ^,x o )e-^- is )l ( _s,-o(7)- 



We remark that W o n cj (— 00, — £) [p(uj, ., .)) G C([0, +00), F) thanks to (16 1. Moreover, p(w, .) G 
C([0,+oc),F) since < +00, where Aj(uj) is defined by ( fTo| and satisfies ( [TI] ). We can 

check that the norm of the operator $ w is bounded by 

||$ w || < if(£,e0 sup sup \V(x,z)\. 

z£[0,L/e] x£[0,d] 



Consequently, using (55), lim £ ^o P(Id — y/e^ is invertible) = 1. Then, the condition (Id — ^fe$ w is 
invertible) is satisfied in the asymptotic e — > 0. On the event (Id — y/e^ is invertible), we have 

W o n w (-oo, -0(p(w, ■)) = (Id - yfe^Y\yfep(u, .) +Po(w, ■)) 

= V^P^, •) + Po(w, .) + T/e$ w (p (v, •)) 



+ OC 



(29) 



5^(v^ w ) J (v^p(w s •) + V~e^(Po(u, ■)))• 
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Moreover, 



||e w o n w (-oo, -0(p{u, ., .)) - v^p(w, .) - po(«, .) - v^* w (po(w, .))IU,f 

< 2e||<i>l ||p( w , .)||oo,f + 2e||$ w || 2 ||h>(w, Olkf 

<if(£,<f)e sup sup |K (x, 2r)| 2 sup \\a(uj,z)\\ n ^ + \\b(w,z)\\n», 

z<£[0,L/e] x£[0,d] z£[0,L/e] 

and therefore 

6 W o iU-oo, -C)(p(w, ., ■)) = Vep(w, ■) + ■) + Ve^(po(w, .)) 

+ C(e sup sup |V(x, z)| 2 sup ||a(w, z)\\ H * + \\b(u, z)\\h? 

V z£[0,L/e] x£[0,d] z£[0,L/e] 

in C([0, +00), FJ . Now, we consider 



P 7 ,2(w,z) 



fc2(w) ^ [ £ ^#(2 ! (^^L/ £ )e i ft" + 6,( W ,zM/e) e -*'') 



( } C ^f (ay(uj, z A L/e)e l ^ u + 6y (w, z A L/e)e-^ u )] d 7 'e^ (,I - 2) du 



fc 2 (w) 



/ E ^^ 3 7P(3i(w 1 zAi/e)e < A« + 6 I ( W> zAL/e) e - < A«) 

JzA,L/e 1 ; = 1 V Pi 



2 V |tI AaV 



k {ay (to, z A L/e)e w ^ u + b Y (u, z A L/e)e^^ u )\ d^e^^du 



Vz e [0,+oo). Using ([T7f, ([T8|, ([T9f, (|20j), and ([29f, we get 

.) -p2{w, -)\\oo,f < K(£,d) y/e sup sup \V(x,z)\ 2 

ze[0,L/e] xe[0,d] 



x( sup \\a{uj,z)\\ H1 > + \\b(u , z)\\ nf + \\p{w,z)\\ F J 



•z£[0,i/e] 

and then 

e w o n w (-oo, -0 (p{uj, ., .)) = ^p 2 (^, .) + po(w, .) + V^0po(^, .)) 



e sup sup |V(x,s)| sup ||a(w,z)|| W |. + + 

z£[0,L/e] xe[0,d] ze[0,L/e] 



inC([0,+oo),.F). Consequently, we can rewrite ( (l7| , ( 18] ), pj| , and (pO) in a closed form in x TC£ . 
Vz € [0,L/e], we get 

^a(w,z) = ^H aa ( W ,z)(S( W ,z))+^H a6 ( W ,z)(i( W ,z)) + Vir is ( W , Z ) 



rfz 



+ eG aa (w,z)(a(w,z)) + e G ab (w, z) z)) + £ R°' s (uj,z) 
+ 0[e 3/2 sup sup |U(>,z)| 2 sup ||a(w, + z)||w" + ||p(w, z)|| F ) , 

V ze[0,L/e] £cS[0,d] zG[0,L/e] ' 

b{u>,z) = ^U ba (uj,z)(a(u Jl z))+^~eU bb (u J ,z)(b(u Jl z))+^~eR b ' Ls (u J7 z) 
+ e G ba (u, z) (a(u, z)) + e G bh (cj, z) (&(w, z)) + eR. ' is (w, z) 

+ o(e 3 / 2 sup sup |U(x,z)| 2 sup \\a(w, z)|| w? + \\b(uj, z)\\ Hf + \\p(uj, z)\\ F ) . 

K z£[0,L/e]xe[0,d] z£[0,L/e] ' 

Let us recall that these equations hold on the event [Id — y/e is invertible) which satisfies the 
condition lim £ ^ W(Id - y/e is invertible) = 1. In these equations, H aa (w,z), H ab (uj,z), H ba (w,z), 
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H bb (uj,z), G aa (uj,z), G ab (uj,z), G ba {uj,z) and G bb (uj,z) are operators from to itself denned by: 



.,2/ \ JV ( W ) C'^iA 



H°>, 2) (y)=H*V> *)(</) = 





JV(u>) 



;=i 

y fe2 (-) C^ 7 ,(*) 



i(A(^)-V^) z 



2 L TTtAM 



Hf ( w> *)( y ) = H>,z)( y ) = ^[ £ .^f we -^("Hft("))» 



H?( W ,*)(y) = H^( J/ ).^[ JV £ 



+ 



fe 2 (w) fY ( \ 

77 ,. ^-v/e~ i(v ^ 7+ ^ )z rf7 / 



7 l/4yV4 y T 



Gf(c;, Z )(i,) = Gf(a;,z)(y) = 



r 



.'7/ 



+ 



ifc 4 (w) 



./£ J -co LjQ 



t^m-VIt'K^-")^ 



+ 



Ve^,(,)C 7V ( U ) 



7"M-VlT'l( u_2 )rfu 



G^(c,z)(y) = G 7 V^)(y) 



Z=l 



+ 



ifc 4 (w) 



|.fe (w) /■-£ z-2 

if J -co *-Jo 



+ 



e -v7 77 «-Vl7'l(2-«)^ u 

/« J -co 



(30) 



(31) 



(32) 



(33) 



(34) 



(35) 
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Gf (w,z)(y) = G* a (u,,z)(y) = 



bar 



L/£ C7 7 ,(«)c^,(«) 



i/3 ! (w)«-^/|7'|(«-z)^ 



cfy e 



yi 



(36) 



+ 



ifc 4 (w) 



fc 2 (w) 



+ 



l7'l(*-«) du 



G°W)(y) = G£V,*)(l/) 



/ 



i = 

C^,{z)C^{u) 



e -*\/7 77 u--y/|7'|(u-z) ( ^, 



d 7 e 



I -i0j(w)z 



- % Pi ( w ) m - x/jyi ( « - u ) ^ u 



y 77' ^^fi^^ -jfl (q,) M - y^vT (u-z) dn 



d 7 'e~^ 



(37) 



ifc 4 (w) 



fc 2 (w) 



Jy/llYlVl* 



L 



i/e ^yWCl »(«) 



+ / " ' ' ' ' _ p-Vt^-VMO'-*),/,, 



The operators H aa (w,2;) and H ab (o>,2;) represent the coupling between the propagating and the ra- 



diating modes with themselves, while the operators Q aa {uj,z) and G ab (w,z) represent the coupling 
between the evanescent modes with the propagating and the radiating modes. Moreover, R a ' is (w, z), 

R a ' s (uj, z), R 6 ' is (w, z), and R ' s (bj,z) represent the influence of the evanescent modes produced 
by the source term on the propagating and the radiating modes. These terms are defined by 



Rf Ls (w,z) = n°^(uj,z) 



b,L s 



lk2{uj) f ^ C ?^<h-,{uj,x )e-^*- L ^die- i ^ 



R"' Ls (^,z) = R 7 < is (w,z) = 



2 ./ , 

ik 2 {u) f-e cy y (z) 

| 7 |l/4 ' 



/: 



(38) 
(39) 



Rf Ls (tj,z)=R:f a (u,z) 



ijfc 4 (w) 



oo */ — 5 



o vWl 



T »(w,a;o)e-VlT"l(«-i S ) e -Vl7'l(^-«) dM 



(40) 



+ 



I/£ C? Y (z)C», Y ,(u) 



7 7 



-^(w,a;o)e-v |7 " l( "- Ls) e-vl 7 'l( u -^rf 



rr 

'-oo J-S 



C»,(z)C», y ,(u) 



L /< C^,(z)C», y ,(u) 

y/VW\ 



„{u,x )e-Vh"\(u-Ls) e -Vh'\(*-u)du 



(41) 
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3.4 Forward Scattering Approximation 

In this section we introduce the forward scattering approximation, which is widely used in the lit- 
erature. In this approximation the coupling between forward- and backward-propagating modes is 
assumed to be negligible compared to the coupling between the forward-propagating modes. We refer 
to [HI E] for justifications on the validity of this approximation. 

The justification of this approximation is as follows. The coupling between a right-going propa- 
gating mode and a left-going propagating mode involves a coefficient of the form 

-|-oo 

E[C£(0)qKz)] cos ((AH + pj(u))z)dz, 





and the coupling between two right-going propagating modes or two left-going propagating modes 
involves a coefficient of the form 

+OC E[quo)qKz)]cos ((AM - 



V(j, i) 6 {1 N(w)} . Therefore, if we assume that 



-t-oo 



E[c;(o)c;(z)] cos ((ftH + ftH)z)rf z = o vc/,0 e {i,...,jvm} 2 , 

then there is no coupling between right-going and left-going propagating modes, which justifies the 
forward scattering approximation, but there is still coupling between right-going propagating modes 
which will be described in Section [4] 



In our context the operator R(uj), introduced in Section 2.1 has a continuous spectrum and it 
becomes technically complex to apply a limit theorem for the rescaled process (a(u>, z/e),b(u>, z/e)). 
The reason is the following. This process is not bounded and the stopping times which are the first 
exit times of closed balls are not lower semicontinuous for the topology of C([0, L],H^ w ), where Wj? w 
stands for TCg equipped with the weak topology. In our context the continuous part (£, k 2 (uj)) of 
the spectrum imposes us to use the norm to control some quantities. Moreover, according to 



Theorem 4.1 in which the energy of the limit process is not conserved, it seems not possible to show 
a limit theorem on C([0,L], (7ijt\ in view of (27 1. In [7j and |9] there is a finite number of 

propagating modes, so that the weak topology and the strong topology are the same. In [11] the 
number of propagating modes increases as e goes to 0. However, in this last case, the problem can 
be corrected by considering the first exit times of a closed ball related to the weak topology and by 
considering the process in an appropriate finite-dimensional dual space. 

In our context if we forget these technical problems, according to [7J |5] the forward scattering 
approximation should be valid in the asymptotic e — > under the assumption that the power spectral 
density of the process V, i.e. the Fourier transform of its z-autocorrelation function, possesses a 
cut-off wavenumber. In other words, we can consider the case where 

nC$(0)C%(z)]coa({ft{u)+l3 j (u))z)dz = V(j,Z) G {l, . . . ,N(ui)} 2 . 

a 

Let us remark that the continuous part (0, k 2 (uj)) of the spectrum, which corresponds to the radiating 
modes, does not play any role in the previous assumption. The reason is that the radiating part of 



the process plays no role in the coupling mechanism as we can sec in Theorems 4.1 and 4.2 below and 
therefore remains constant. 

Finally, we shall consider the simplified equation on [0, L/e], 

a{u),z) = ^H aa {cj,z){a{u,z))+^K a ' Ls {uJ,z) 



dz 



+ e G aa (uj, z) (o(w, z)) + e lT' Ls (w, z) 



+ o(e 3/2 sup sup \V(x,z)\ 2 sup \\a(u,z)\\ n « + \\p(u, z)\\ F ) 

K »6[0,i/e] xe[0,<J\ ze[0,L/e] 7 

in Ti^ . We shall see in Section |4j under the forward scattering approximation, that 

lim P (L e < L) = VL>0, 
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where 



U 



inf(i>Q, sup ||3(w,2!)||^ + ||p(w,z)|||.>4=) 

ze[0,L/e] ~ ^/ £/ 



Consequently, we can show that V 77 > 



limP( sup ||a(u;,z)||L. - ||a (w)||^ > 77) = 0. 



sup 

' z£[0,L/e] 



This result means that the local energy flux for the propagating and the radiating modes is conserved 
in the asymptotic e — > 0. 



4 Coupled Mode Processes 

In this section, we study the asymptotic behavior, as e — > in first and £ — * in second, of the 
statistical properties of the coupling mechanism in terms of a diffusion process. 
Let us define the rescaled process 



a £ (w, z) — a yw, -J Vze[0,L]. 



This scaling corresponds to the size of the random section [0, L/e]. This process satisfies the rescaled 
coupled mode equations on [0, L] 



+ G" («, J) (a>, *)) + R Q ^ S (a;, f ) (42) 
+ 0(v / e sup sup |1/(a;,z)| 2 sup \\a e (uj, z)\\ n ^ + \\p(uj, z/e)\\ F ) 

K z€[0,L/e] x£[0,d] ze[0,L] ' 

in Tig, with the initial condition a e (ui,0) — So(w). We shall see that under the forward scattering 
approximation the condition L e > L is readily fulfilled in the asymptotic e goes to 0. 

Proposition 4.1 VL > 0, 

limP(L e < L) = 0, 

w/iere 

L 

«£[0,i/e] 



inf(i>0, sup \\a(uj,z)\\2 +\\p(uj,z)\\ 2 F > ^-) : 



lim lim P sup ||a £ (w,z)||i„ > M = 0. 
This result means that the amplitude a e (w, z) is asymptotically uniformly bounded in the limit e — > 



on [0, L]. More precisely, according to Section 3.2 we have V 77 > 



lim . 



'( sup \\o?{uj,z)\\^ - ||ab(«)||^« >v)=0, 



that is the local energy flux for the propagating and the radiating modes is conserved in the asymptotic 
e -> 0. 

Proof Using Gronwall's inequality, VL > we get 

lim Tim" P ( sup \\a?(u, z)||L, > M, L < L e ) =0. 

M^+ooe^O \ze[Q,L] 5 / 
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This result means that the process a E (uj, .) is asymptotically uniformly bounded on [0, L] and then L e 
is large compared to L in the asymptotic e — > 0. In fact, VX > and VM > 



' (i e < L) < P L e < L, sup ||a e (w, z) ||^ < M 

\ 2G[0,LAL 5 ] 5 , 



Moreover, 



+ P sup \\a e (LU,z)\\^ > M 
\ze[0,LAL>] 5 . 



L e < L, sup ||a e (w,z)||!^ < M = 



se[0,LAi £ 



for e small enough, since for L e < L 



s- 1/2 < sup +||p( W ,.)|| 
ze[o,L e ] 5 



<M + K(e,d)e sup sup \V(x,z)\ 2 M + 2\\p (< J j,.)f OOjF 

zS[0,L/e] ze[0,d] 



according to (29 1 



Let us introduce a\ (w, .) the unique solution of the differential equation on [0, L] 
-a\(uj,z) = i=H QQ (u, (at(u,z)) + G aa («, J) (a!(w,«)) 



<iz 



(43) 



in 7i^, with initial condition o,\(lj, 0) = ao(o;). Using Gronwall's inequality and (54l we can state that 



lim lim 



sup \\a\(uj,z)\\n- > M = 0. 
v ze[o,L] / 



The relation between the solution of the full system ( 42 ) and the one of the simplified system ( 43 1 is 
given by the following proposition. 

Proposition 4.2 



V?7 > and V/i > 0, limP sup ||a e (cj, z) - a\(uj, z)\\ H » > r\ = 0. 

^° \zeln.L] 5 / 



Proposition 4.2 means that the information about the evanescent part of the source profile is lost in the 



asymptotic e goes to 0. In fact, the coupling mechanism described by the system (42 1 implies that the 



information about the evanescent part of the source profile is transmitted to the propagating modes 



through the coefficients R a ' Ls (uj,z) and R""~ a (w,z) defined by (|38f, (j39l), Bfjf and Ell. In these 



- a.L 



expressions we have the term 7 /(o;,a;)e vIt'K 2 is ) which comes from the right-hand side of ( 24 1 



and which is the form of evanescent modes without a random perturbation. This term is responsible 
for the loss of information about the evanescent part of the source profile because of its exponentially 
decreasing behavior. 

Proof We begin by proving that VL > 0, V77 > and Vp > 



limP sup \\a e (LU,z) - a\{oj , z)\\ 2 H ^ > 77, L < L e = 0. 



In fact, R a ' Ls (u>, z) decreases exponentially fast wit h the propagation distance. Moreover, R a (u>, z) 
can be treated as G aa in the proof of Theorem 4.1 be caus e e _ V / M('"~ i s) canno t be compensated by 
e -«/3j(")z nor by e - l V^ z . Moreover, using Proposition 4.1 we get the result. ■ 
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Finally, we introduce the transfer operator T^' e (o;, z) from TCg to itself, which is the unique operator 
solution of the differential equation 

^T«< e ( w , z) = ±=n aa (u, J) T«->, z) + G aa (u, ^) T*' e (w, z) (44) 
with T e ' e (w,0) = Id. Then, 

Vze [0,£], oi(w,z) = T«' £ (w,z)(a (w)), 
and we get the following result. 
Proposition 4.3 

Vr? > and > 0, limPI sup ||a e (w, z) - T^' e (w, z)(a (w))||^ > ?/ | = 0. 
4.1 Limit Theorem 

This section presents the basic theoretical results of this paper. In [9] and [14], the authors used 
the limit theorem stated in |22j since the number of propagating modes was fixed. However, in 
our configuration, in addition to the _/V(w)-discrete propagating modes the wave field consists of a 
continuum of radiating modes. The two following results are based on a diffusion-approximation 
result for the solution of an ordinary differential equation with random coefficients. This result is an 
extension of that stated in [32] to the case of processes with values in a Hilbert space. 

Theorem 4.1 VX > andVy efi^ C^xtf^Pfu)), the family (T ? ' e (w, .)0/)) ee (o i)> solution 
of the differential equation ( 44 1 , converges in distribution on C([0, L],W ffl ) as e — * to a limit denoted 



by T 4 (w, .)(y). Here H£ w stands for the Hilbert space equipped with the weak topology. This limit 
is the unique diffusion process on Tit, starting from y, associated to the infinitesimal generator 



ru ru I rid , ru 



where 



N(uj) 



£ i = \ E r M w ) (T j T j &r l djr + T l T l &r j d T T-T j T l &r j d Tl - ^1}^%) 

3,1=1 

i^H / _ _ x 

+ 2 E r ^ w ) (Wt^ + T^^ - TjTid Tj d Tl - TjVdjrdjr) 

3,1=1 

2 e cw - r ») (' A + 5?%) + i e w - 55%) . 



2 

3=1 3=1 



and 



c lz 2 E ( A f H + lA f M) / '-' / ' + ( A f M - ' A 7 ; ^'0 

3=1 

AT(w) 



3=1 



Here, we have considered the classical complex derivative with the following notation: If v = v\ +iv2, 
then d v = \ (d Vl — id V2 ) and dp = h (d Vl +id V2 ). We have used the following notations. V(j, /) S 
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{1,..., JVH} 2 andj^/ 



/ \ /-+00 

r M^) = oflv W? s / E[qi(o)q{(*)] cos (09, (w) - &H»<fe, 



2pi{w)h{u>) Jo 



2=1 



1,4/, ,\ r+oo 



W = o^.tl ^ / E[qi(0)C5{(z)] sin ((AH - ^(w))z)dz, 



2/3»AM J 



and V(i,0 S {l,...,iVH} 2 , 



1.4/ \ /+oo 



2pi{w)ft{u) Jo 

/•fe 2 (") 1.4 /-A /+oo 

AfH = jf 2Sm(dJ E[C^(0)C^(*)]ooB((VV-/3i(a;))*)d«*/, 
2^/(0,) 7 E[^(0)^(z)]sin(( v ^-/5 J H)z)dzd 7 ', 

^ 1,4/- \ /+oo 

k*Au>)= y / E[C^,(0)C; y (z)]cos(/3,Hz)e-Vl7'l^^ 7 '. 



The coupling coefficients C w (z) are defined by (15 1. We get the following result in the asymptotic 

Theorem 4.2 VL > ond Vy € 7Yq = x L 2 (0, fc 2 (») ; tfie /ami/?/ (T«(w, .)(z/)) ee(0 1} converges 

in distribution on C([0,L], (TCq, ll-llwg')) as £ — * to a limit denoted by T (u>, .)(?/). This limit is the 
unique diffusion process on Hq, starting from y, associated to the infinitesimal generator 



3 ■ 



whe 



£ 2 = ~\ E ( A i( w ) + iA ») % + ( A » - iA j( w )) T i%> 

3=1 



3=1 

Here, we have Vj £ { 1, . . . , N(u) } 

A5(w) = ]imA? c (w), A?(w) = Hm A*' e (cA <«.•(«) = lim/cffw). 

Theorems |4.1| and |4.2| describe the asymptotic behavior, as e — > first and £ — > second, of 
the statistical properties of the transfer operator T^' c (o;,L), in terms of a diffusion process. In the 
appendix we give the proofs of Theorems |4.l| and |4.2| which are based on a martingale approach using 
the perturbed-test-function method. In a first step we show the tightness of the processes, and in a 
second step we characterize all subsequence limits by mean of a well-posed martingale problem in a 
Hilbert space. 

The infinitesimal generator C u is composed of three parts which represent different behaviors on 
the diffusion process. We can remark that the infinitesimal generator depends only on the N(ui)- 
discrete coordinates. The first operator describes the coupling between the _/V(w)-propagating 



C. GOMEZ 



SHALLOW-WATER PROPAGATION 21 



modes. This part is of the form of the infinitesimal generator obtained in [7J|9], and the total energy 
is conserved. The second operator L% describes the coupling between the propagating modes with 
the radiating modes. This part implies a mode-dependent and frequency-dependent attenuation on 



the iV(w)-propagating modes that we study in Section 5.1 and a mode-dependent and frequency- 
dependent phase modulation. The third operator £3 describes the coupling between the propagating 
and the evanescent modes, and implies a mode-dependent and frequency-dependent phase modulation. 
The purely imaginary part of the operator does not remove energy from the propagating modes 
but gives an effective dispersion. 

Moreover, let us remark that the convergence in Theorem 4.1 holds also on C([0, L], (Tt^ 
for the A r (w)-discrete propagating mode amplitudes. 



4.2 Mean Mode Amplitudes 

In this section we study the asymptotic mean mode amplitudes. From Theorem |4.2| we get the 
following result about the mean mode amplitudes. 

Proposition 4.4 Vy E H% , Vz E [0,L], Vj € {l, . . . , N(u)} 



lim lim E 



Tf e (cj,z)(y)] =E[T$(w,z)(y) 



exp 



Z + I 



y 3 {w). 



(45) 



First, let us remark that the mean amplitude of the radiating part remains constant on £ 2 (0, k 2 (uj)). 
Second, Vj S {l, . . . , N(iv)}, the coefficient (T^j(uj) + Aj(lu) — T^Auj))/2 is nonnegative. In fact, 
for 6 {1, . . . , N(oj)} 2 such that j ^ I, Tj^ui) and rj^(w) are nonnegative because they are 

proportional to the power spectral density of and at (3i — 0j (u>) and frequencies. Therefore, 
— Tjj(u>) is also nonnegative. Moreover, A^(w) is also nonnegative because it is proportional to the 
integral over (0, k 2 (u;)) of the power spectral density of at ^7 — (3j(uj) frequency. 
The exponential decay rate for the mean jth-propagating mode is given by 



T>,L)(y) 



y,- exp 



r^.M-r^H + A^ 



j j 



L 



which depends on the effective coupling between the propagating modes, and the coupling between 
the propagating and the radiating modes. This exponential decay corresponds to a loss of coherence 
of the transmitted field. 



5 Coupled Power Equations 

This section is devoted to the analysis of the asymptotic mean mode powers of the propagating modes. 
More precisely, we study the asymptotic effects of the coupling between the propagating modes with 
the radiating modes. Let 



Tjtuj.z) = lim limE 

J C^O e^O 



|Tf(u,,Z)(</)| 



E 



|T?(^)G/)| 



(46) 



be the asymptotic mean mode power of the jth propagating modes, ^(uj, L) is the expected power 
of the jth propagating mode at the propagation distance z = L. Here y l E Hq is defined by = Sji 
and yl^ — for 7 E (0, fc 2 (w)), and where 5ji is the Kronecker symbol. The initial condition y l means 
that an impulse equal to one charges only the Zth propagating mode. From Theorem |4. 2 1 we have the 
coupled power equations: 



-Aj(u)5>,2) 



N(uj) 

71=1 



(47) 



C. GOMEZ 



SHALLOW-WATER PROPAGATION 22 









z 






. ' Ocean 




d \ 






Bottom 



Figure 3: Illustration of the radiative loss in the shallow-water random waveguide model. 



with initial conditions TJ(uj,0) = Sji. These equations describe the transfer of energy between the 
propagating modes and T c (lo) is the energy transport matrix. In our context, we also have the 
coefficients Aj(u) given by the coupling between the propagating modes with the radiating modes. 
These coefficients, defined in Theorem |4. 2 [ are responsible for the radiative loss of energy in the ocean 
bottom (see Figure [3]) . This loss of energy is described more precisely in the following section. 



5.1 Exponential Decay of the Propagating Modes Energy 

In this section, we assume that at least one of the coefficients A c (uj) is positive. With this assumption, 
we show that the total energy carried by the propagating modes decays exponentially with the size L 
of the random section. In the opposite situation, that is when there is no radiative loss A c (lj) = 0, it 
has been shown in [7] and [HJ Chapter 20] that the energy of the propagating modes is conserved and 
for large L the asymptotic distribution of the energy becomes uniform over the propagating modes. 
Let us defined 



X j > Vj G {1, .. . , N(lu)} and \\X\\* M = (X,X). 



1 



with (X,Y) t 



ES x j Y j for {X,Y) G (I"'"*) 2 , and 

A§(w) = diag(A c 1 (uj),...,A c N{u) (uj)). 



Theorem 5.1 Let us assume that the energy transport matrix T c (tu) is irreducible. Then, we have 

1 



lim — In 

L — >+oo L 



i=i 



= -Aoo(w) 



with 



AooH = inf ((-r» + A^H)x,x) 



xes 



(48) 



+ 



which is positive as soon as one of the coefficients A c Auj) is positive. 



This result means that the total energy carried by the expected powers of the propagating modes 
decays exponentially with the propagation distance, and the decay rate can be expressed in terms of 
a variational formula over a finite-dimensional space. 

Proof The coupled power equations admit a probabilistic representation in terms of a jump Markov 
process. If we denote by (l^ ) t>0 a jump Markov process with state space {1, . . . ,N(lu)} and 
intensity matrix r c (w), then we have using the Feynman-Kac formula: 



Thcj,z) = E 



exp 



A'' 



)(w)(is I 1 



= I 



(49) 
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Moreover, we have supposed that T c (lu) is irreducible. Then, (lj in an ergodic process with 

invariant measure (J-n(u)i which is the uniform distribution over {1, . . . , N(u>)}. That is, fJ-N(u>)(j) = 
1/N(u>) Vj € {1, . . . , N(w)}. The self-adjoint generator of the jump Markov process (Yj^ )t>o is 
given by 

N(u>) 
n=l 

for every function (/> from {1, . . . , N(uj)} to M, and it is easy to check that £ N ^ /j>n(u) — 0- Let us 
consider the local times 



o 1 K M =i) 



f/.S 



for j e {1, . . . ,N(u>)} and T > 0, which corresponds to the time spent by the process (Y t ) t>0 in 
the state j during the time interval [0, T]. According to [5], we have a large deviation principle for 
^It viewed as a random process with values in M± which is the set of probability measures on 
{1, . . . , N(oj)}. More precisely, we have 



lim -InE 

L^+oo L 



exp 



lim — In E 

£— H-oo L 



exp 



(- f A c Y N ( ^ds^ 



- inf (7( M ) + (A C M, A1 )) 



with 

Consequently, 



i(fx) = || ( - r c H) 1/2 ^||2 RtfM = {(- r») Vm, Vm), 



lim — In 



JV(u>) 



-Aoo(o;). 



Let us assume that A oc (cl') = 0. As S^^' is a compact space, there exists X e 5^"'' such that 

AocCw) = (( - r» + A5(w))X , a ) rnm = 0. 

Moreover, —T c (lu) and AS(w) are two nonnegative matrices and is a simple eigenvalue of — r c (w) by 
the Perron-Frobenius theorem. Then, 

((-T c (uj))Xo, X ) RN[LO) = X = y/fiff( u ), 

and 

(AS(w)Xo, A ) RJVM = => 3j e {1, . . . , N(w)}, X (j) = 0. 

Therefore, 

Aoo(uj) > 0. 



The expression (48) of Aoo(w) is not simple. However, we have the following inequalities. 



je{i,...,N{u)} 



A' 



'M < A 00 (w) < A(w) = £ A?(w) 

.7=1 



(50) 



First, we assume that Vj € {1, ... , AT (a;)}, A|(o;) = A(a>) > 0. In this case, using (50) 

Aoo(a;) = A(w). 

This means that if all the coefficients which represent the radiation losses are equal, the decay rate 
of the total energy of the propagating modes is given by this coefficient. 
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Second, we assume that the coupling matrix is small, that is, we replace r c (w) by tT c (lu) with 
t < 1. If Vj e {1, . . . , N(oj)}, A^lu) > we have 



H A^(w) = 



je{i,...,iV(w)} 



A5(«). 



From (50 1, it is the smallest value that Aoo(o;) can take. This result is consistent with the fact that 
the coupling process on the transfer of energy between propagating modes is negligible and the decay 
rate of the energy of a particular propagating mode j is given by its own decay coefficient Aj(u>). 
Then, for the total energy of propagating modes the decay rate is given by the minimum of those 
decay coefficients. Consequently, if there exists A c ^ q (lo) = 0, we have 

limA^H = 0. 

r— +0 

The reason is the energy of the joth propagating mode stays approximately constant with a weak 
transfer of energy, and 

1 



lim -A^H = inf ((-r»)X,X) RNM > 0, 



where 



t^o r xev 

N(u 



V = {xe suppX c {1, . . .,N(oj)} \ su PP (A c (uj))} , 



V. 



because ^JHn(u., 

Now, we assume that the coupling matrix is large, that is we replace T c (oj) by ^T c (oj) with t<1. 
In this case, we have 

lim A^(w) = A(w). 



From (50 1, it is the largest value that A oc (u;) can take. The strong coupling produces a uniform 
distribution of energy over the propagating modes and the decay rate becomes (A c (w), [J,n{w)/™n(u) = 
A(w) for each mode. A more convenient way to get this result is to use a probabilistic representation. 
In fact, we have 



Tj(u,z)=E 



E 



exp 



exp 



A! 



Af 



Y, 



N(ui) 



= I 



where [Y, 



>t>0 



is a jump Markov process with state space {1, . . . ,N(u)} and intensity matrix 



r c (w). Using the ergodic properties of (Y t ) t>0 , we get that 



lim T T,l (ui, L) 



N(oj) 



exp (-A(w)i) 



Finally, if we assume that the radiation losses are negligible, that is we replace A c (w) by tA c (u>) with 
r < 1, we have 

limA;M = 0. 

In fact, if the radiative loss is negligible, the coupling process becomes dominant, and we can show 
that 

) = 0{t), 



VL>0, sup \\TJ^z)-T^{uj,z)\\ %] 

ze[Q,L] 



where T°' l (ui, .) satisfies (47 1 without the coefficient A c (oj). In this situation 

T°' l (u,L) = P A 



N(u) 







and the total energy is conserved (see Figure |4| , and 

lim-A^H = AH>0. 

r^O T 

As it was already observed in [5] the modal energy distribution converges as L 
distribution: 

lim 7?'\u,L) = — --. 

£-»+oo J N{U!) 



-oo to a uniform 
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Figure 4: Illustration of negligible radiation losses in the shallow-water random waveguide model. 



5.2 High- Frequency Approximation to Coupled Power Equations 



In this section, under the assumption that nearest-neighbor coupling, introduced in Section 5.2.1 



is 



the main power transfer mechanism, we give an approximate solution of the coupled power equations 



(47 1 in the high-frequency regime or in the limit of large number of propagating modes N(u) 3> 1. 



Let us note that the limit of a large number of propagating modes N(u) 3> 1 corresponds to the 
high-frequency regime lo — > +oo. Next, we analyze the energy carried by the propagating modes in 
this regime. 

The coupled power equations can be approximated in the high-frequency regime by a diffusion 
equation. This approximation has been already obtained in |14| for instance, in which we can find fur- 
ther references about this topic. We can also refer to |18j for more discussions on this approximation. 
For an application of such a diffusion model to acoustic propagation in randomly sound channels we 
refer to [T!5] , and for applications to time reversal of waves we refer to [TT] . 

Using the form of the covariance function d52| , we find 



and 



where 



ak 4 (uj)Ij : i(Lj) 



2/3 i (w)AM(o 2 + (/3 J (w)-AM) a ) 



ash 



fc2(w) afc 4 (o;)J 3 , 7 H 

2/?,MV7(a 2 + (/3 J H-V7) 2 ) 



drf, 



Ijl =-A?A? S((Tj - ai,(Tj - CT;) + S((Tj + ai,(Tj + er ; ) 

— S(aj — <tj, (7j ; + (ii) — S((jj + ai , <ta — <tA 



T. —- A 2 A 2 



S ((Tj — T), <Jj — rf) + S(aj + rj, aj + rj) 
- s ( a 3 ~ V, <* i +v)-S (a j + rj, Oj ~ 7 l) 



with 



S(vi,v 2 ) = / / 7o(a;i,X2)cos (-r^i) cos {-^-X2)dxidx2 
Jo Jo ad 



and where Aj(u>), A 7 (u>), <Jj(w), r](u>), (j)j(u),x), and 7 (w,x) are defined in Section 2.1 



5.2.1 Band-Limiting Idealization 

In this section, we introduce a band-limiting idealization hypothesis in which the power spectral 
density of the random fluctuations is assumed to be limited in both the transverse and the longitudinal 
directions. 
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We assume that the support of S lies in the square [ — ] X [ — , ] ■ Then, 

J i otherwise, 



and 



— ) i j 



l o otherwise. 

From this assumption we get VO < 7 < k 2 (w) and j G {!,..., N(ui) — 2}, 



rj(uj) — (Jj(u>) > nik(id)di 1 ^ — <T j( Ct; ) > n\k(uj)d9 — (N(u>) — 2)n 

>n( nik ^ d e-N( U ))+2n. 



e[o,i) 



Then, for j E {!,... , N(u>) - 2}, 



3tt 



inf — > 



and 



A 3 c H = o, v. 7 e{i,...,7VH-2}. 



Consequently, the coupled power equations (47 1 become 
^(z) = -A%Tl(z) + T%_ 1N (T^iz) - T l N (z)) 



"t~^/V-i( 2: ) — 7 ^A r -l'^/V-l( z ) + T^-l N-2 (^JV-2( Z ) — ^/V-l( z )) 



+ T< k _ 1N (T<,(z)~T^_ 1 (z)) 1 (51) 

- r^-u Ctf-i W - ^W) + r i + i,' - for i e {2, . . . ,iV - 2}, 

^T/(z)=r= 1 (T 2 i (^)-T 1 i (^)), 
with 7/(0) = fy. 

The band-limiting idealization hypothesis is tantamount to a nearest-neighbor coupling. More 
precisely, this assumption implies that V(j, I) € {1, ... , N(uj)} 2 the jth mode amplitude can exchange 
informations with the Ith amplitude mode if they are direct neighbors, that is, if they satisfy \j — l\ < 1. 

5.2.2 High-Frequency Approximation 

The evolution of the mean mode powers of the propagating modes can be described, in the high- 
frequency regime or in the limit of a large number of propagating modes N(u>) 3> 1, by a diffusion 
model. This diffusive continuous model is equipped with boundary conditions which take into account 
the effect of the radiating modes at the bottom and the free surface of the waveguide (see Figure [3] 
page^ 



221. 



Let, G C°([0,1]), Vw G [0,1], and z > 0, 

j=l 

where (p 1— ► T^^\z,.) can be extended to an operator from L 2 (0,1) to itself. Here, L 2 (0, 1) is 
equipped with the inner product defined as follows: V(ip,ip) G L 2 (0, l) 2 

(V'V^rafoji = / (p(v)ip(v)dv. 
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Theorem 5.2 We have 

1. Mip G L 2 (0,1) andVz > ; 

Urn T"W(z,u)=T v (z,u) m£ 2 (0,l), 
where T v (z,u) satisfies the partial differential equation : V(z,u) G (0,+oo) x (0,1), 



wi£/i the boundary conditions 

d 

—77,(2,0) = 0, T^U, 1) = 0, ond 77,(0, u) = <p(u), 
ou 

Vz > 0. 

2. Vw € [0, 1], Vz > 0, ond \/tp G C°([0, 1]) smc/i that tp(l) = 0, we have 

lim Tj v ^(«,u) = 7;(« s «). 

w— ►+00 ^ 

-Here, 

aoo(u) 



l-(l-^)(^) 2 ' 

with ao — 2a ^ifd°e 2 > ®~ ~~ V n i - ^0 = Jjf 70(^1: ^2) cos (f^i) cos {^2X-f)dx\dx2- n\ is the index 
of refraction in the ocean section [0,d], X/a is the correlation length of the random inhomogeneities 
in the longitudinal direction, and 70 is the covariance function of the random inhomogeneities in the 
transverse direction. 

This theorem is a continuum approximation in the limit of a large number of propagating modes 
N(lo) 3> 1. This approximation gives us, in the high-frequency regime, a diffusion model for the 
transfer of energy between the N{uj) -discrete propagating modes, with a reflecting boundary condition 
at x = (the top of the waveguide in Figure [2] page [3| and an absorbing boundary condition at u = 1 
(the bottom of the waveguide in Figure [2| which represents the radiative loss (see Figure [3]) . 

5.2.3 Exponential Decay in the High-Frequency Regime 

In the high-frequency regime, we also observe that the energy carried by the continuum of propagating 
modes decays exponentially with the propagation distance. The exponential decay of the energy in 
the high-frequency regime is given by the following result. 

- 2 



Theorem 5.3 Vip G L 2 (0, 1) \ {0} such that ip>0, and Vu G [0, 1), 

lim i In [77, (Z, it)] = — A^, 

L— >+oo Li 



whe 



1 

2, 



and 



Aqo = inf / a^v)^'^) dv > 
Jo 

T> = U>ec°°([o,i]), |MU 3( o,i) = 1, !^(o) = o, v(i)=o 

This result means that the energy carried by each propagating modes decays exponentially with the 
propagation propagation, and the decay rate can be expressed in terms of a variational formula. 
Consequently, the spatial inhomogeneities of the medium and the geometry of the shallow-water 
waveguide lead us to an exponential decay phenomenon caused by the radiative loss into the ocean 
bottom. 
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Proof We can see that the operator P^ = g^(floo(')g^) on ^ 2 ([0j 1]); with domain 

V(Poo) = |ye H 2 (0,1), A^(0)=0, ¥>(!)= 

is self-adjoint. has a compact resolvent i?> = (Xld — Poc)^ 1 because [0, 1] is a compact set and 
then it has a point spectrum (Xj)j>i with eigenvectors denoted by {<fioa,j)j>i- Moreover, all the 
eigenspaces are finite-dimensional subspaces of T>{P oa ) and Vp £ V^P^) \ {0} 

Let us organize the point spectrum in the nonincreasing way, • • • < A 2 < Ai < 0. We have 



%p( l ,v) = ^(v3,0 OO j) L 2 (Oil) e A3i 0o O j(w). 
j">i 

Lemma 5.1 Ai is a simple eigenvalue and one can choose </>oo,i such that </>oo.i( u ) > Vt> £ [0, 1). 



Proof (of Lemma 5.1) This lemma is a consequence of the Krein-Rutman theorem, but not its 
strongest form |2B]. Indeed, the set of nonnegative functions in L 2 ([0,1]) has an empty interior. 
However, using the smoothness of the eigenvector the proof also works in our case as we shall see it. 

Using the maximum principle we know that if <p> £ L 2 ([0, 1]) such that ip > 0, we have %p(L, .) > 0, 
and then R\(ip) > 0. Consequently, applying the Krein-Rutman theorem [26 to the resolvent operator 
R\ with A > and which is a compact operator, the spectral radius p(R\) is an eigenvalue, for which 
one can associate an eigenvector p p (r x ) such that Vw £ [0,1], P p (r x )(v) > 0. However, we have 
Vt> £ [0, 1), <pp(R x ) (v) > 0. In fact, let us assume that there exists vq £ [0, 1) such that P p (r x ){vq) = 0j 
then Rx(ipR x )(v ) = 0. Moreover, R x (pr x ) = p(R\)tp p ( Rx) is an eigenvector for and then P p {r x ) 



is a smooth function on [0, 1]. Therefore, according to the proof of Theorem 5.2 we have 

f+OO 



I) 



R\{<PRx)(vo) = I e x % HRx) (t,v )dt 

[p p{Rx) (x(t))l {t<Tl) ]dt 



g-AtT, 







1 e- Xt p p(Rx) (\x(t)\)dt 



0. 



where P„ is the unique solution of the martingale problem associated to = (5oo(")J^) anc ^ 
starting from vq. Here, we have chosen such that Vw £ [0,1], aoo(v) — a oc (— v) — aoo(v), and 
the martingale problem associated to £ Soo is well-posed. Moreover, t\ = inf(i > 0, \x(t)\ > 1). 
Consequently, P„ ( JJ" 1 e~ M [ ip P (R x )(\x(t)\)dt = 0) = 1. However, we know that there exists V\ £ (0, 1) 
such that <pp(R x )(yx) > 0, and then v% < vq < 1. Therefore, P„ (ti < t Vi ) — 1, and by the Markov 
property 

< E p "o [e- T »a (Tvi<+oo) ] = E p "o [e- r "il (Tui<+00iTl<Tui) ] 

<E Pl [e- T "H (T „ i<+oo) ] 
<E F "o [e-»H K<+co) ], 

which is impossible. Therefore, £ [0,1), P p (r x ) > 0. Now, to see that the eigenvalue p(R\) is 
simple, let <p £ L 2 (0, 1) \ {0} such that R\(tp) — p(R\)ip, and let 

P Rx :R — > C°([0,1]) 
i i — ► - tip, 

which is a continuous function. We recall that ip is a smooth function on [0,1]. Let us show that 
3t £ K such that 9? = tip p / Rx ^, that is € Pr a (R). To do this let us assume that £ Pr a (R). By 
linearity one can assume that 3fo € [0, 1) such that p(vq) > 0. Let 77 > be small enough to have 
v £ [0,1 - rf\. Let Jf+ = {ip £ C°([0,1 - r/]),V« € [0, 1 - 7?],y>(u) > 0}, then the interior of K+ 
for the sup norm on [0, 1] is K++ = {p £ C°([0, 1 - r/]),Vu £ [0,1 — rj\,<p(v) > 0}. Moreover, for t 



C. GOMEZ 



SHALLOW-WATER PROPAGATION 29 



small enough p Rx — tp> E A'+ + , and p Rx — tp £ for t large enough. Then 3t E K such that 
— toV G K£ \ However, p Rx — top > 0, but pn x — t ip ^ because ^ Pr x (R). Following 

the previous work we have 

p(R^(p Rx - t p>) = R\(<p Rx - t p) E K+ + . 
Consequently, p{R\) = 1/(A — Ai) implies that Ai is also a simple eigenvalue and one can choose 

0oo,l = R\(fRx) = P{R\)V>Rx S K+ + . 

That concludes the proof of Lemma |5.1| □ 

As a result, Mp S L 2 (0, 1) \ {0} such that p > 0, Vw E [0, 1) we get 

lim y In [T V (L, v)] — Ai, 

and 

Ai = sup (PooW), <p) L 3/ [0 m = ~ A oo < 0. 

( ( 3GX'(Poo) 



IIVllz,3([o,i]) 



=1 



In Theorem 5.3 we take cp E L 2 (0, 1) \ {0} such that cp > 0, which can be consider as being the initial 
repartition of energy over the continuum of modes. However, the result of Theorem |5.3| is also valid 
for any p E L 2 (0, 1) \ {0} such that (<p, <?W) L 2 (CM) > 0. 

5.3 High-Frequency Approximation to Coupled Power Equation with Neg- 
ligible Radiation Losses 

In the case of negligible radiation losses, we also get a continuous diffusive model for the coupled 
power equations in the high-frequency regime or in the limit of a large number of propagating modes 
N(lu) 1. This diffusive continuous model is equipped with boundary conditions which take into 
account the negligible effect of the radiation losses at the bottom and the free surface of the waveguide 
(see Figure [3] page 25 ) . 



Now, let us assume that the radiation losses are negligible, that is, A c (lu) = tA c (lu) with t« 1. 
We have already remarked that, if the radiation losses are negligible, then the coupling process is 
predominant and we have 

VL>0, sup ||77V,2)-^°'W)lk R "("> =0(t), 
ze[o,L] 



where T°' 1 (oj, .) satisfies 



d 

dz 

±-T?\z) = V Uj - T°- l (z)) + T° j+lj (7&(*) - Tf'izj) for j E {2, . . . , N 1}, 

dz 1 



with 7^(0) = 8ji. 

5.3.1 High Frequency Approximation 

Let, Vp> E C°([0, 1]), Vu S [0, 1], and z > 0, 



where p i— > T^ u \z, .) can be extended into an operator from L 2 (0, 1) to itself. 
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Theorem 5.4 We have 

1. Vip G £ 2 (0,1) andVz > 0, 

lim T*^(z,u) = T tp (z,u) inL 2 {Q,\), 

where T v (z,u) satisfies the partial differential equation : V(z,u) G (0,+oo) x (0,1), 

d d ( d 

with the boundary conditions 

d_ 

du VK ' " ; dv 

Vz > 0. 

2. Vm G [0, 1), Vz > 0, and Vip G C°([0, 1]) such that <p(l) = 0, we have 



-T v (z,0) = 0, — T v {z, 1) = 0, and T v (0,u) = tp(u), 



Here, 



lim T"M{z,u) = T v {z,u). 



aoo(w) 



1 - (1 - (^) 2 

ittttfj a = 2a ^4g2 , = V 1 - l / n i, So = J d /o d 'T ( Xl ' a;2 ) cos ( 7 3 Xl ) cos (f x 2 )dxidx 2 - n x is the index 
of refraction in the ocean section [0,d], l/a is the correlation length of the random inhomogeneities 
in the longitudinal direction, and 70 is the covariance function of the random inhomogeneities in the 
transverse direction. 

This theorem is a continuum approximation in the limit of a large number of propagating modes in 
the case where the radiation losses are negligible. This approximation gives us, in the high-frequency 
regime, a diffusion model for the transfer of energy between the 7V(tj)-discrete propagating modes, 
with two reflecting boundary conditions at u = (the top of the waveguide in Figure [2] page [3]) and 
u = 1 (the bottom of the waveguide in Figure |2j. Here, the two reflecting boundary conditions mean 
that there is no radiative loss anymore (see Figure [4]) . 

5.3.2 Asymptotic behavior of T(L, v) as L — > +00 

In the case where the radiation losses are negligible, we have seen in Section [O] that the decay rate 
satisfies lim r ^o ^L( tJ ) = an d T ' l (w, L) converge to the uniform distribution over {1, . . . , N(ui)} as 
L — > +00 [5]. In the high-frequency regime we have the following continuous version. 

Theorem 5.5 V<p G L 2 (0, 1) and Vtt G [0, 1], 

lim T v (L,u) — / ip(v)dv, 

that is, the energy carried by the continuum of propagating modes converges exponentially fast to the 
uniform distribution over [0, 1] as L — > +00. 

As a result, the energy is conserved and the modal energy distribution converges to a uniform distri- 
bution as L — > +00. 

Proof We can see that the operator P^ = ^(^oo(O^) 011 ^ 2 (p7 1]); with domain 

^gf 2 (o,i), q v <p(o) = o, ^^( 1 ) = ° 

is self-adjoint. Moreover, P^ has a compact resolvant because [0, 1] is a compact set and then it has 
a point spectrum (\j)j>a with eigenvectors denoted by (<f>oo,j)j>o- Moreover, all the eigenspaces are 
finite-dimensional subspaces of r D(P OQ ) and V<£ G r D(P 00 ) \ {0} 

( P °°(^)>¥>}l2 (0 ,i) 
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Let us remark that Ao = is a simple eigenvalue with eigenvector </>oo.o = 1- Then, the spectrum is 
include in (— oo,0] and we have the following decomposition 

T<p{z,v) = I ^(w)c?t; + ^(( y 5 ! OO j) i2(Oil) e A ^o3 oo , J (t;). 
i>i 



o 



Therefore, Vit € [0, 1], 

with exponential rate Ai < 0. 



lim T v (L,u) — I tp(v)dv, 

L— >+oo 



Conclusion 

In this paper we have analyzed the propagation of waves in a shallow- water acoustic waveguide with 
random perturbations. In such a waveguide, the wave field can be decomposed into three kinds 
of modes, which are the propagating, the radiating, and the evanescent modes, and the random 
perturbations produce a coupling between these modes. 

We have shown that the evolution of the propagating mode amplitudes can be described as a 



diffusion process (Theorems 4.1 and 4.2 1. This diffusion takes into account the main coupling mecha- 
nisms: The coupling with the evanescent modes induces a mode-dependent and frequency-dependent 
phase modulation on the propagating modes, the coupling with the radiating modes, in addition 
to a mode-dependent and frequency-dependent phase modulation, induces a mode-dependent and 
frequency-dependent attenuation on the propagating modes. In other words, the propagating modes 
lose energy in the form of radiation into the bottom of the waveguide and their total energy decays 
exponentially with the propagation distance. We can express the decay rate in terms of a variational 



formula over a finite-dimensional space (Theorem 5.1 ) 



Under the assumption that nearest-neighbor coupling is the main power transfer mechanism, the 
evolution of the mean mode powers of the propagating modes can be described, in the high frequency 
regime or in the limit of a large number of propagating modes, by a continuous diffusive model with 
boundary conditions which take into account the effect of the radiation losses at the bottom and the 
free surface of the waveguide. In this regime, we observe that the energy carried by the continuum of 
propagating modes also decay exponentially with the propagation distance. The exponential decay 



rate can be express in terms of a variational formula (Theorem 5.3). 

The diffusive systems obtained in this paper can be used to analyze pulse propagation and refo- 
cusing during time-reversal experiments in underwater acoustics [TT] . 



6 Appendix 

6.1 Gaussian Random Field 

This section is a short remainder about some properties of Gaussian random fields that we use in the 
proofs of Theorems |4.1( and in Sections |3.2| and |3.3| All the results exposed in this section can be 
shown using the standard properties of Gaussian random fields presented in pQ and [5] for instance. 

In this paper, the random perturbations of the medium parameters are modeled using a random 
process denoted by (V(x,t),x € [0, d],t > 0). Throughout this paper the process V is a continuous 
real-valued zero-mean Gaussian field with a covariance function given by 

E[V(x, Zl )V(y,z 2 )} = l0 ( x ,y)e~ a ^-^ V(x,y) e [0,d] 2 and W(z u z 2 ) S [0,+oo) 2 . (52) 

Here, a > 0; 70 : [0, d] x [0, d] — > R is a Lipschitz function, which is the kernel of a nonnegative operator, 
that is, there exists a nonnegative operator Q 7o from L 2 (0, d) to itself such that V(</5, ip) <E L 2 (0, d) 2 

r-d r-d 

(Qto(f),^) L 2 i0 _ d) = J J jo(x,y)(p(x)ip(y)dxdy. 

Consequently, one can consider the process (V(.,t))t>o as being a continuous zero-mean Gaussian 
field with values in L 2 (0, d) and covariance operator Q l0 . In other words, Vn G N*, V(<£>i, . . . , (p n ) € 
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L 2 (0, d) n , and V(t a . . . , t n ) e [0, +00)™ 

(V^(*i), ■ • ■ , V Vn (t n )) = ((V(., ti),Vi) L2{0>d y. ■ ■ , (V(.,t n ), ^„> L2(0)d) ) 
is a real-valued zero-mean Gaussian vector such that V(j, /) S {1, . . . , n} 2 

E[V^fe)^(*0] = (Qio^ 3 ),w) L , [04 f- a]t] - tl] - (53) 
With this point of view we have the following proposition. 
Proposition 6.1 We have 

1. (V(.,t)) t >Q is a continuous zero-mean stationary Gaussian field with values in L 2 (0,d) and 



autocorrelation function given by (53). Then, we have Vn £ N* and Vi > 0, 

I'd \ 



E 







\V(x, t)| dx 



= E 







|F(x,0)| dx 



< +00. 



(54) 



2. We have the following Markov property. Let 

F t = a{V(.,s),s<t) 
be the a-algebra generated by (V(., s),s < t). We have 

(V(.,t+h)\r t ) = (v(.,t + h)\a(V(.,t))), 
where the equality holds in law, and this law is the one of a Gaussian field with mean 

E[V(.,t + h)\T t ] =e- ah V(.,t) 
and covariance, \/((p,ip) € L 2 (0,d) 2 , 



E 



-2ah\ 



V v (t + h) (t + h) - E [V v (t + h) I ft] E [V^ (t + h) \ T t ] 

= (Q lo ( v ),^) L2(0d) (1-e 

The Markov property of the random process (V(.,t))t>o is a direct consequence of the exponential 
form of the autocorrelation function (53 ) with respect to the variable t [1]. This property will be used 
in the proof of Theorems |4.1| which are based on the perturbed-test-function method. 

Now, we are interested in some estimation on the suprcmum of V(x,t) with respect to the two 
variables x and t. To this end, let us introduce some notations [2J. Let e > be a small parameter 
and L > 0. We consider the following pseudo-metric on the square [Q,d] x [0,L/e] defined by 

m((x, t), (y, s)) = E [(V(x, t) - V(y, s)) 21 1 2 



< K l0 [\t 



y\}- 



Let us remark that [0, d] x [0, L/e] associated to the pseudo-metric m is a compact set. From Theorem 
1.3.3 in 0, we have 



E 



sup \V(x, t) 

te[0,L/e] 



< K 



diam([0,d]x[0,L/e])/2 



H 1/2 {r)dr 



™p ie [0,d] 7o(x,x) 



/In K 



dL 



dr. 



where H(r) = ln(iV(r)), and N(r) denotes the smallest number of balls, for the pseudo-metric m, 
with radius r to cover the square [0,d] x [0,L/e\. Here, diam stands for the diameter with respect to 
the pseudo-metric m. Consequently, we have the following proposition. 

Proposition 6.2 V/x > and WK > 0, 



limPfe^ 1 sup sup \V(x,t)\>K 

^ x£[0,d] te[0,L/e] 



0. 



(55) 



Moreover, according to Theorem 2.1.1 in |2J, one can show that the limit (55 1 is obtained exponentially 
fast as e — > 0. 
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6.2 Proof of Theorem I3~T1 

The proof of this theorem is in two parts. The process (T^' £ (z)) z>0 is not adapted with respect to 
the filtration T% — !F z / e . Then, the first part of the proof consists in simplifying the problem and 
introducing a new process for which the martingale approach can be used. The first part of the proof 
follows the ideas of [15]. The second part of proof of this theorem is based on a martingale approach 
using the perturbed-test-function method and follows the ideas developed in 0]. 

Then, let us introduce T (.) the unique solution of the differential equation 

^-T Le (z) = -^H aa (-) ¥'\z) + (G aa )T 4,e (z), (56) 
dz v e V e / x ' 

with T ? ' e (0) = Id and where (G QQ ) is defined, Vy € W s , by 

/-£ ,-1,4 f+oo 
— — / E[C jy (0)C, V (*)] cos (prie-VM'dzdr/yj 
-oo 2Pj^\j'\ Jo 

Vj € {l, . . . , iV} and (G aa )^(y) — for 7 <E (£, k 2 ). We have the following proposition that describes 
the relation between the two processes T^ £ (z) and T^' (z). 
Proposition 6.3 

Vy €Ht andVr]>0, limP ( sup \\T^(z)(y) - (z^y)^ > r, ) = 0. 

\ze[o,L] J 

Let us remark that the new process (T^' (z)) >q is adapted to the filtration !F* and 

\\T te (z)(y)\\ 2 Hi = \\yf Ui Vz>0. 

Let r y = \\y\\n ( , 

Br y ,n ( = {a e \\\\\ Ht = ^{\\) Hi < r v} 

the closed ball with radius r y , and {g n -,n > 1} a dense subset of Br y ,*H*- We equip B ry- n i with the 
distance dB TyiH defined by 

V(A,/i) € (B ry ,-H £ ) 2 , and then (Bji i ,dis ry H ) is a compact metric space. 

Using a particular tightness criterion, we prove the tightness of the family (T^' Q)ee(o,i) on 
C([0, +00), (B rvi -H(: d-B ry:7i )), which is a polish space. We have chosen such a space to be able to 
apply the Portmanteau theorem. In a second part, we shall characterize all subsequence limits as 
solutions of a well-posed martingale problem in the Hilbert space TL^ . 

We have the following version of the Arzela-Ascoli theorem [T3] for processes with values in a 
complete separable metric space. 

Theorem 6.1 A set B C C([0, +00), (B ry _-H ( , dts ry H )) has a compact closure if and only if 

VT > 0, lim sup tot(<7, rj) = 0, 

»7^° g£A 

with 

m T {g,T])= sup d B (g(s),g(t)). 

0,t)6[0,T] 2 

\t-s\<n 

From this result, we obtain the classical tightness criterion. 
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Theorem 6.2 A 

and only if 



family of probability measure (P e ) eg( - ^ on C([0,+oo), (B. 



r v ,Ht, a B r 



) ) is tight if 



VT>0,?/>0 lim sup ¥ e (g ; m T (g, v) > v') = 0- 



From the definition of the metric d& ry H , the tightness criterion becomes the following. 

Theorem 6.3 A family of processes (^ e ) c e(o,i) is tight in C([0, +oo), (B ry ,u s , de ry n )) if and only 
if ((X\ X) n ) ee{0 1} is tight on C([0, +oo),C) VA G ^. 

This last theorem looks like the tightness criterion of Mitoma and Fouque [?T1 fH] . 

the family 



For any A G 7i^, we set T^' £ (z)(?/) = (T q '"(z)(y), A) w ■ According to Theorem 



6.3 



(T 4 ' C (.)(y)) £ is tight in C([0,+oo), (B ry! n v d Bry ^)) if and only if the family (f^' £ (.)(y)) £ is tight on 
C([0,+oo),C) VA G H ( . Furthermore, (T (.)(y)) e is a family of continuous processes. Then, it is 

sufficient to prove that VA G Ti.^, (T A ' (.)(y)) e is tight in 2?([0, +oo), C), which is the set of cad-lag 
functions with values in C equipped with the Skorokhod topology. 



Proof (of Proposition 6.3) Differentiating the square norm and using the fact that H aa (z) is skew 
Hcrmitian, we get 



< 2 



\T^(z)(y)-f"\z)(y)f Ht 



2\\(G aa )\\ [ Z \\T^(u)(y)-T^(u)(y)\\ 2 ni . 
Jo 



Let rf > 0, we will split the interval [0, z/e] into intervals of length rf j \fe. The idea is that over these 
intervals the fast dynamic of G aa averages out while T^' e does not move significantly. We have 

: f Z/ \{G aa (u)) - (G aa ))T^%eu)(y),T^(eu)(y) -¥-\eu)(y)) du 
Jo v ' n i 



< 



G aa {u) (G aa ))T^{eu){y),T^{eu){y) - T ^ {eu){y)\ du 

I Hi; 

f l ((G°» - (G aa ))T^(eu)(y),T^eu)(y)-f £ '\eu)(y)) n du 



with 



7?J s/i 



f G aa (u) - (G aa ))T^(eu)(y),T^(eu)(y) - T 5 'V)(y)\ du 

I Tic 



< 



e 1/4 v/V 



G c 







(") 



1/2 



'du 



vWH(G aa )| 



x sup \\T^(z)(y)\\ n A\T^(z)(y)-¥'\z)( y )\\^ 
«e[o,£] 



since < z — [yf^y] v^ 7 / < v^ 7 /; an d 



> G aa {u) „ ( G -»T^(eu)(t/),T«' e (e U )(j/) - T f 'V)(y)\ du 



V 1_I M m + 
m=0 Jm Vi 



(m+l)4= 



;G»»(«) - (G aa ))T 4,c (eu)(y), T^' e (ew)(y) - f ? ' £ ( ra )( 2/ )\ d U 

' He 
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Moreover, 



and 



(eu)(y) = T^{mr,'sfe){y) + / VeH aa (v)T^(ev)(y) + eG aa (v)T^(ev)(y)dv 
{eu){y) = ¥'\mr)'Ve){y)+ /" , ^H aa ( W )T e ' e (6 V )(y) + e(G aa )T^(e W )(y)dt;. 



Therefore, we have 

,(m+l)^ 



^_ ,(m+l) ^ / ^ _ (G aa ))T ^ (eu)(y))T «,e M(2/) _ tf>* (eu)(y)\ du 



V '7/ ((G a » - (G aa ))T«' e (mr ? '^)(y),T«' e (mr ? '^)(y) - f ^(m V 'V~e)(y))du 



/ , / , e((G aa (u)-(G aa ))H aa (v)T^(ev)(y),T^(eu)(y)-T ( ' (eu)(y) 

+ e 3 / 2 ((G aa ( U )-(G aa ))G aa ( U )T«- e M( 2 ;),T^(6 U )(y)-T^M(y) 



+ e((G aa ( U )-(G aa ))T^(mr 7 'yi)(y),H aa ( W )(T^(e V )(y)-T ? ' e M( 2/ )) 



+ e 3 / 2 ((G a » - (G aa ))T^ £ (mr7'^)(y),G aa ( W )T^ £ (e W )(y) - (G aa )T C ' e ^)(y) ) du. 



Consequently, by the Gronwall's inequality 



sup \\T^(z)(y) - T-(z)(y)\& ( < B(e, V ')e 2 
ze[o,L] 



\(g-)\ 



where 

B(e,rf) = 2 



x sup \\T^(z)( y )\\ n( \\T^(z)(y)-T-(z)(y)\\ ni 

ze[o,L] 



r^ri-i 



+ 2V~e £ / , / , (24l|G aa H|| + ||(G-)ll]H Haa M 



m=0 



+ e 3 / 2 [||G a >)|| + ||<G Qa )||] 2 ) sup \\T^%z)(y)\\ n( \\T^%z)(y)^T^(z)(y)\\ n( dvdu 

' ze[o,L] 



+ 



'' I {(G aa («) - (G aa ))T«< e (mr/^)(y),T^(W^)(2/) - T^(mr/yi)(y) > du 



and 



( sup ||T^(z)(j/)-t €,e (z)(y)||^ >„) < p(s(e, t,') > ^- 2 IK G ")H L ). 
y ze[o,L] 



Setting rj" = -qe 



Kg-) | 



we have 



•(B(e,»/) > r/') < P(s(e,r ? ') > ry", sup \\T^(z)(y)f H( < m) 
+ P( sup ||T^(*)(y)|& 4 >M). 
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We already know that the process T ' (.)(y) is bounded. Moreover 



B(e,r)') > V ", sup ||T^ e (z)(y)||^ < M) < — E 



B(e, V ')l 



(™P*S[0,.L] \\^''(z)(y)\\ 2 n <m) 



with 



E 



B(e,rj')l 



(sup.elo.iillT 4 ''^)^)!!^ < M ) 



< K[if + e 1 /^ + Vifa' + r/ 2 )] 



[^]- 



2^ ^ IE i 1 (sup i!e[0 , I , ] |tT«.«(^)( y )||^ <M) 



m— 
(m+1)-^ 



<X[r ? ' 2 + e 1 /4 v ^7 + ^ + ^2 ) ] 
+ 2^feK E 



(G aa (u) - (G oa »T«' £ (m )? '^)(y),T^ e (m» ? 'Vi)(y) - T e 'W^)(y)) ^ 



m=0 



(m+l)^= 



G aa (u) - (G aa )du 



1/2 



since f / ^ ^ f" n' dv du — — and 



E 



G° 



(~) f] <KE (f \V(x,0)\ 2 dx) 



(57) 



for u € [0, L]. As a result, it remains us to estimate only one term. 
Lemma 6.1 



lim Je V E 



m— 



2 

G aa (u) - (G QQ )dw 



1/2 



Proof (of Lemma 6.1) Let us remark that we have the following decomposition. For each j S 
{l, . . . , iV}, almost every 7 6 (£, fc 2 ), and My € Wf, 

G° a (z)(y) = ^G^ a (z)y ; + / G^,(z)yyd'y', 
1=1 

G™(z)(y) =Y,G°?(z) yi + / G~,(*)t^dY. 

7 1 C 



Letting 



/■(m+l)5= 

/ G aa (u) - (G aa )d Ul 

V e 
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wc have (j, I) 2 e {l, . . . , iV} 2 such that j ^ I, and almost every 7 <E (£, k 2 ) 



P = 

17 33 



Jl 



P / 

r 31 



P , 



















= Ve 


ImK. 






= Ve 










,(m+l)^ 







G%(u)-(G™)..du, 

G°°, (u)du, 
G a J(u)du, 



G a ",(u)du, 



and 



JV JV fc 2 fe 2 N ,2 

\wn 2 <Ei p ^i 2 +E / Ei p ^i 2 rf7+/ / 



1 2 j„/ 



c?7 c?7. 



Moreover, 



E 



F(xi,zi)t/(x2,Z2)V r (x3,z 3 )y(a;4,Z4) =£ V(ari, Zi)V(a:2, z 2 ) E V(x 3 , z 3 )V(xi, z 4 ) 

+ £ [V(si , «i) Vfo, z 3 )j £ [v(x 2 , 22)^(2:4, z 4 ) 
+ E [V(si , zi) V(ai4, z 4 )] -E [v(a; 2 , za) V(ar 3 , 23) 
- 7o(*i, £2)70(^3, a: 4 ) e — I- 1 — le—l^-^l 



-a\zi-Zi\ e -a\z2-z 3 \ 



+ 7o(zi, ^4)70(^2, x 3 )e 



which is the fourth order moment of a Gaussian field. To compute the expectation of the square 
norm of P we must know these moments. Following that decomposition, the square norm of P can be 
decomposed in three parts. First, after a long computation, the two parts corresponding to the two 
last terms of the previous decomposition are dominated by y/e uniformly in m. Then, we focus our 
attention on the part corresponding to the first part of the previous decomposition. For E[|Pj 7 /| 2 ] , 
E[|P 7 ;| 2 ] , and E[|Pj;| 2 ] with j ^ I, we get after a long computation terms of the form 



/ / e i ^-^ Ul e- i ^-^ U2 du 1 du 2 = 0(e), 

,(m+l)4 ,(m+l)< 

e/ / e l(/3i ~V7)«i e -^(A-V7)«2 duidu2 = 

/ ^ / ^ eM-M^e-M-M^d^ du 2 = 0(e). 



(. — "p= "lit: f— 

rKl)^ /•(™+l)^r 
f / / 

For E[|P 77 /| 2 ] we separate the integral into two parts 



f f E[|P 7 y| 2 ](i7'd7= f E[|P 7 y| 2 ]d 7 '(i7+ f E[|P 7 y| 2 ]d 7 'd7, 



where \x > and 



^ = {(7,7')e(£,fc 2 ) 2 , |V7-V7|>m}, 
^ = {(7,7')e(C,fc 2 ) 2 , \Vi-VY\<v}- 
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Consequently, 



f E[|P 7 y| 2 ]d 7 'd7 < K f dj'dj, 



and on J> M we get terms of the form 



(m+l)^= /•(m+l)3 



AVy-^i e -*(VY-^2 dui du 2 dj'dj = 0(e). 



Now, it remains us to study E[|Pjj | 2 ] . After a long computation, the terms of order one produced 
by G a ^ are compensated by the terms of order one given by (G aa ) .. Moreover, the other terms are 
dominated by ^/e. 
As a result, we get 



lim J~t V E 



m=Q 



(m+1)^ 

y/l I , G aa (u) - (G aa ) du 



1/2 



< if J / d 1 'd 1 



and one can conclude the proof of Lemma [6. 1| by letting /j, —> O.D 
From the previous lemma, we finally get, V?/ > 



lim P( sup \\T^(z)(y)-T 6 ' £ (z)(y)f Hi > 



1(^)1 



V) < 



-Kri 



12 



since using the Gronwall's inequality and ( 57 1 we have 



lim lim p( sup ||T 4 ' e (z)(j/)||^ > M ) = Q 



26[0,L] 



Consequently, we conclude the proof of Proposition 6.3 by letting 7/ — > 0. 



According to Proposition 6.3 to study the convergence in distribution of the process (T^' c (.)(?/)) e 

it suffices to study the convergence for (T ' (•)(#))• Moreover, we shall consider the complex case 
for more convenient manipulations. Letting A € TL^, we consider the equation 

with H x = (H,X) n , G\ = ((G aa )(.),X) H , where, for j € {l, . . . , N} and almost every 7 e (£, fc 2 ) 



fc 2 /-~Y 



2 TTtA 7 i/4yV4 

The proof of Theorem 4.1 is based on the perturbed-test-function approach. Using the notion of 
a pseudogenerator, we prove tightness and characterize all subsequence limits. 



6.2.1 Pseudogenerator 

We recall the techniques developed by Kurtz and Kushner [IB]. Let .M e be the set of all J rc -measurable 
functions f(t) for which sup t<T E [\f(t) \] < +00 and where T > is fixed. The p — lim and the 
pseudogenerator are defined as follows. Let / and f 5 in A4 e V<5 > 0. We say that / — p — lima f 5 if 

supE[|/ 5 (£)|] < +00 and lim E[|/ 5 (£) - f(t)\] = Vt. 
t,s 
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The domain of A e is denoted by T>(A e ). We say that / 6 T> (A e ) and A e f = g if / and g are in 
V(A e ) and 

" Ef [/(* + *)]-/(*) gW 



n — lim 

5^0 



0. 



where Ej is the conditional expectation given ,F t e and J 7 ^ — T t / t . A useful result about ^4 6 is given 
by the following theorem. 

Theorem 6.4 Let f e V(A e ). Then 



is an (J^t) -martingale. 



M}(t) = f(t) - / A e f(u)du 



6.2.2 Tightness 

We consider the classical complex derivative with the following notation: If v — a + i(3, then d v = 
\ (d a - idp) and dy — \ (d a + idp). 

Proposition 6.4 VA € Ti^, the family (T A (•)(i')) ee (o i) * s tyht on 2? ([0, +oo),C). 

Proof According to Theorem 4 in [IB], we need to show the three following lemmas. Let A £ H^, f 
be a smooth function, and /o(t) = / (t a (t)(y)j . We have, 



+ chf (T^(t)(y) 



1 
1 

7e 



t\ t 



H X [T- (t)(v),C[- l,-)+G A (T s ' (t)(y) 



F A (T ? ' £ (t)(y),C 



t\ t 
ej'e 



Ga fT ? ' £ (t)(y) 



Let 



f!(t) = — a„/ (t1' £ W(z/)) jf E t e [h a (f c ' e (t)(y),c (H) , 

+ (f f(i)(y)) y +0 °E t e [ff A (f C,e (t)(y),C (") , ") 



du. 



Lemma 6.2 VT > 0, lim c sup 0<t<T |/i(i)| = almost surely, and sup t>0 E [|/f(i)|] = O (\/e). 
Proof (of Lemma 6.2 ) Using the Markov property of the Gaussian field V, we have 



with 



a,/ Tr (t)(») Ff,(t) 



^.aW = E 



E 



1 n\t)(y) 



a + i{(3i - fa) 
« 2 + (Pi ~ PjY 



4^ 



« 2 + (^-/3,) 2 



A, 



V- grj (j) e i(ft-^y)j T^fflM a + »(A-V7) 



c 



+ (Vy - x/7) : 



X-ydj. 
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Using (54 1, we easily get 
and 



E[|/i £ (i)|] <V~eK(f,\). 



\f!(t)\<K(X,f)Ve sup sup \V(x,t)\ 

0<t<T/e xG[0,d] 



Then, we can conclude with (55).D 

Lemma 6.3 {A € (/„ + ft) (*), e € (0, 1), < t < T} is uniformly integrable VT > 0. 



Proof (of Lemma 6.3) A computation gives us 

A<(ti + ft)(t) = F x (T 6 ' e (t)(y),C r ' 



where 



for (qi,q 2 ,q 3 , <? 4 ) € ({1, . . . , TV} U (£, k 2 )) , with 



C 



C(T) £g> C(T) qi q2 q3 qi — C qi q2 (T)C q3q4 (T) 
4 



t \ t 



and 



F x (T,C,s) = d v f(T) F^(T,C, S ) + G X (T) + ^/(T) i^T, C,s) + G x (T) 



d 2 J(T)F 2 (T, C, s) + 9|/(T)F A 2 (T, C, s) 



(kd v f(T)F 3 x (T, C, s) + 9^/(T)F A 3 (T, C, s), 



3=1 



in.ptin.. V a 2 + _ p.yi 



c 



^ a 2 + (A _ /3 . ) 2 d T 



7 ^^-(Vt 7 -^) 2 



fc_ 4 /■'• 

4 



N r k 2 
,k 2 N 



7 2 



« 2 + (A - Vt) 



a + i(VY-\/l) ^ 
a 2 + (VY- V^)' 



4 5 VVry'vr 



C 77Y^ ei (^^)s T a + - ft) 

a +(V7'-\/7) 2 



A 7 c?7, 
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^a 2 (T,C, S ) 

= 4E 

lil ^ a 2 + (A-/3 J ) 2 72 



V Cjij'l' -iQg,-i8j+j9,,-j9 J ,) J r r ^ a + *(0l ~ Pj) 

^NPWPv a2 + (Pi ~ Pj) 2 



+ 



E 

1=1 



'JTiJ'*' P »(\/7i-ft+ft'-/?j/)8 T , T 



« 2 + (y7[-/5,) 2 " 



d 7 i 



4 W fe 2 



"7i 7 2 „ 2 



a 2 + (VVi-ft) S 



W ,-fc 2 



E 



c* 



/=1 ^ \Jp~jP^V%l2 
,fc 2 N 

/ s = 



72 a 2 + (A - fy) 



<*7i 



1 2 a 2 + (v / 7T-ft') 2 



\jp3^fiil2ll 



k 



4 /-fc 2 JV 



x jrt x JTnPiPWv i a 2 + (A _^ r) 2 



w ,fe 2 



y VriPiPj' \fi 2 



: « 2 + (A-VtT) ; 



,fc 2 JV 



^ i'=i yjyFtHJyiv 71 a 



a + *(V7i - v/Ti) , , 



+ (VtT- VtT) 



fc 2 ,fc 2 a 



y/Sttb>SZ 72 « 2 + (v^-V7r) 2 7172 



A 7l Aj/ 



JV 



E 



Z=l 

,fe 2 JV 



/* C 7i7(72<' p »(^-V7T+ft/-V7i)« T ., T;( i 

^ i'=i Jy/irfsnfa 71 a 



d 7l 



+ (vVi - a/tT) 



7« A (7W2 7 2 ) 1/4 71 72 a 2 + (y 7 [-v / 7T) 2 



A 7l A 72 rf 7 id 7 2, 
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N 

E 



V Cjlj ' 1 ' p i(.Pi-Pi-Pi>+M s T,TZ a + ^ - ^ 

a + i(/3 ( - /3,) , 

"72 



72 a2 + (A-^)- 



r fc 2 # r _ 



l v a /' f= r"' 

71 a*+( V^-ft) 2 



r k2 , 2 C 



71 72 « 2 + (v/7[-^) 2 



A, A 



•j A?' 



E 

7=1 



E 



c 



: jh2i' cl (f3 l -f} 1 -p„+^)s T] 7jr; « + j(A - A;) 

L/i ^a 2 + (A-/3i) 2 ' 2 



i 'a 2 + (A-/3 J ) 2 
d 7 ' 



^ — 71 a 2 + (V7i-/3 J -) 2 ' 



7 ? ('=1 



— a + i(V7i -A) , , 



71 72 a 2 + (v^[-/3,) 2 



AjA 72 (i72 



__ 
4 



E 



' + (A - yfn) 



•E, 



r k 2 N r 

/ 7i7 1 j''' e i(V7r-V7r-ft/+/V)»T -TV — 

if ^ . / /^7«..fl.. 71 " 2 



1.2 .1-2 



i - l72 a 2 + ( / 3 i -v7r) 2 " 72 

7^7i 



« 2 + (771-771) 



71 72 a 2 + (v^-V7T) 2 



A 7l Aj' 



.4 /-fc 2 /-fc 2 



AT 



q + »(A - VtT) 
2 + (A - VTi) 2 



1 72 a 2 + (A-V7T) 2 72 



I y-> ^7i7i72^ c »( A /7T-^-ft/+^)s T ^.=r; 



— a + ^(77T-V7i) , , 



7 i ' .2 



d 7 l 



« 2 + (v / 7[-y7T) : 



^ /"* C 7l7i727^ _ e ^(^-^7T-^+^)«T^TV ", " ^ _ ^-^d 7 W72 



(7i7W27 2 ) 1/4 



Ti 7 2 2 



A 7l A 72 d7ic?72. 



This expression combined with (54 1 gives us, sup e t E[ |^l £ (/q + f{) (t)\ 2 ] < +oo. □ 
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Lemma 6.4 



lim limP ( sup |fi' e (i)(y)| > M ) = 0. 

M^+oo e^O \ <t<T 



Proof (of Lemma 6.4) We recall that ||T ' (t)(y)||H s = \\y\\H i and then 

|Ti' e (i)(y)| < \\¥' e (t)(y)\\ H( \\\\\ H( = |M| Wt ||A|| Wt . 

□ 

This last lemma completes the proof of Proposition |6.4| ■ 
6.2.3 Martingale problem 

In this section, we shall characterize all subsequence limits by showing they are solution of a well- 
posed martingale problem. To do that, we consider a converging subsequence of (T^' C (.)(y)) c6 ( 01 ) 

which converges to a limit T^(.)(y). For the sake of simplicity we denote by (T (.)(z/)) e e(o,i) the 
subsequence. 

Convergence Result 

Proposition 6.5 VA G H.^ and V/ smooth test function, 
f(T{(t)(y)) 

d v f(T{(s)(y)) (j«(T e ( S )(y)),A) +ckf(T{(s)(y))(j^(s)(y)),x) 



+d 2 v f(T{(s)(y)) (K(T*(a)(y)){\),\) + %f{T*(a)(y)){K{T*(8)(v))(\),\) 



+ckd v f(T{(s)(y)) (L(T«( s )(y))(A),A) + 9«9jt/(t|(s) (y))(L(T^(s)(y))(X),X) ds 



is a martingale, where 



H 



rh -a. / (n, -a* 

+ i 



-J ' 



i=i i=i 

N N 

L(T)(Xh = ,E r j' T J T|Al + o E r i' T ' T ' A i- 



/=i 
'A; 



and 



J«(T) 7 = X(T)(A) 7 = L(T)(A) 7 = 
/or almost every 7 <E (£, fc 2 ), and /or (T, A) € 7i|. 



Proof (of Proposition 6.5) Let 

f +00 



w) 



w, 



F } 



(T ( '\t)(y), 



\(T te (t)(y),E[C(0)®C(0)], 



c{-) •"(")• 



Lemma 6.5 



supE[|/ 2 e (i)|]=0(e) 

i>0 
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.4'(./u -!-/! -KAK/.i: Fx (t €,e (t)(j/),E[C(0)®C(0)]^ ) f.l !<./>. 



wftere sup t>0 E[|A(e,t)|] =0(y/e). 



Proof (of Lemma 6.5) Using a change of variable we get /J(£) = eB(e,t) with 



/ \ (T S,e (*)(y),c( U + t\®c(u+ -) ,u+ f - 



o 

# A (t 6,e (t){y),E[C(Q) <g> C(Q)],u + J ) ,1,,. 



By a computation, we get that sup e t>0 E [\B(e, t)\] < +oo, and after a long but straightforward 
computation we get the second part of the lemma. □ 

Next, let G A (T 5 ' £ (i)(2/), f) = F A (T C ' £ (t)(y),E[C(0) <g> G(0)], f) and 



./:(/)= - /' \G x (T te (t)(y),~) - T ^ go ~ [ T G x (f^(t)(y),s)ds 







du. 



Lemma 6.6 VT" > 0. we have 



Urn sup E[|/|(i)|]=0. 



'0<t<T' 



Proof (of Lemma 6.6) Using a change of variable, we get 



Let /x > 0, we have 



du. 



where for j S {1,2,3,4}, 

= |(7i)ie{i,...,3} € (C, fc2 ) : ' 1 3((7 i )iG{i,...,4-j} G {/3i,---,Av} 



1 <M 



4-j 



3 4-j 

and (w)ie{i,...,4} G {-1, l} 4 , with ^m^fi +^2^ l+] qi < fj,j. 



i=i 



1=1 



Finally, 



lira sup E 

e ^°o<t<r' 





t 






e 


s: 







1 /*■* 

G A (^(ii.uj-^^ ^ G A (T«- £ (i)(y), S i r/.s- 



rfi/ 



4 /■ 

/ d 7l ...d 7j , 

n 1 * Is*,, 



3=1 " J <^ 

and then by letting /i — > we get the announced result. □ 
Let / e (i) = / Q e (i) + fl(t) + /|(f) + /|(i). A computation gives 

i r T 

A c f'(t) = lim - G x (T e>e (t)(y),a)d8 + C(e,t) 
= Gr (T C '\t)(y))+C(e,t), 
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where, V/i > 0, 



lim sup E[\C(e,t)\] <K(T',£,y) V / . . . d 7j , 

3=1 <f 

using the boundness condition ( [54] ). Moreover, for (T, A) € W|, G°° is defined as follow 



Gf(T) =d v f (T) (j«(T), A> w< + ckf (T) (J« (T), \) 



dtf (T) (X (T) (A), A) + d±f (T) (if (T) (A), A 

\ / 7~L£ \ / rip 



chjd v f (T) ( L (T) (A), A) + 3 v (kf (T) (L (T) (A), A 



where 

^(T)(A) 3 - 

L(T)(A), 

for j e {1,...,N}, with 
and 



fc4 



E 



\ - C j7j'F T T a + - ft) 



k 1 ^ C iW rp ffr- ° + *(A - ft) x 

= J- / J- /' — s ^T7?Ai' 



c = E[a(o)®c(o)] > 

if(T)(A) 7 = L(T)(A) 7 = 



for almost every 7 S (£, A; 2 ). 

We assume that the following nondcgcncracy condition holds. The wavenumbers ft are distinct 
along with their sums and differences. This assumption is also considered in [7], [9] and [14] ■ As a 
result we get 



G? (f 6 ' € (t)(y)) = d v f (T^(t)(y)) (jHT^(t)(y)),\ 



chf (Tf(t)(y)) (j€(T e,6 (t)(y)),A 
^/(tl' e (t)(y)) (if (T*' e (i)(y)) (A), A 



a|/ (fi' e (t)(y)) (if (T^(t)(y)) (A), A 
^/(fi' e (t)(j/)) (i(f^(t)(y)) (A), A 



(58) 



He 



He 



d v chf (ri:\t)(y)) (L(T^(t)(y)) (A), A 



He 



By Theorem 6.4 (ilf| e (t)) >Q is an (J 7 /) -martingale. Then, for every bounded continuous function h, 
every sequence < si < ■ ■ • < s n < s < t, and every family (Xj)je{i,...,n} with values in H.^ we have 



E 



h (t%( 8j M, 1 < j < n) (f(t) - f(s) - f A*r(u)d- 



= 0. 



Finally, using (58) and Lemmas 6.2 6.5 and 6.6 we can conclude the proof of Proposition 6.5 



Uniqueness In order to prove uniqueness, we decompose T^(.)(y) into real and imaginary parts. 
Then, let us consider the new process 



Y«(i) 



Y 2 ' ? (i) 



where Y 1,e (i) = fle(T e (i)(») and Y 2 ^(t) = im(T ? (t)(y)) . 
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This new process takes its values in x Q^, where = M. N x £ 2 ((£, k 2 ), M). x is equipped 
with the inner product defined by 



< T > s > e5 x e , = E T } s } + T ' s2 + / 



V(T, S) € <?f X <? s . We also use the notation Y{(t) = (Y c (i), A) with Ae^x We introduce 
the operator T on x given by 



-T 1 



By Proposition |6. 5 1 we get the following result. 
Proposition 6.6 VA S £ 5 x V/ G C, 



6°( 



/(Y|(t))- / <B«(Y<W),A> /(YlW) 



+ i<^(Y<(a))(A) 1 A> ( , ex0 /'(YiW)d B 



is a martingale, where 

with, forj e{l,...,N} 
B Z(Y) 3 = 



A(Y)(A) - A X (Y)(A) + A 2 (Y)(A) + A 3 (Y)(A), 



A' 



2 



Ti(Y) 



A 1 (Y)(A) 3 = T,.(Y) + T 2 (Y)A 2 ] 

i=i 

N N 

A 2 (Y)(X), = —Yj J2 n [YfAj + Y 2 A 2 ] + T, (Y) £ r?, [T, 1 (Y)A i 1 + T 2 (Y)A 2 ] 



z=i 



i=i 



N 



A 3 (Y)(A) J ; = A,]T [(Y, 1 ) 2 + (Y 2 ) 2 ] , 



and 



B« (Y) = A 7 (Y)(A) = A 7 (Y)(A) = ^(Y)(A) = 
/or almost every 7 € (£, k 2 ), and for (Y, A) <E (p£ x ^) 2 - 
Proof (of Proposition |6.6[ ) Using Proposition |6.5| 

/(Yi(t))-^i?e((j«(T«( s )( 2/ )),A) H5 )/'(Y|( s )) 

+ \Re(({L + K)(T* (s)(y))(X), A)^)/" {Y{(s))ds 
is a martingale. Let us remark that we also denote by A the function A 1 + iA 2 , and 

and Jm((T«(t)(y),A)^) = (T(Y«(t)),A) 



Re((Tt(t)(y),\) ni ) = (Yt(t),\) 
Then, we have 

Re((jt(Tt(s)(y)),\) n( ) = <i^(Y«( S )), A>^ 
Re(((L + K)(^(s)(y))(X),X) ) = (A(Y«(,s))(A), A) 



e 5 xe 5 ' 
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As a consequence of Proposition 6.6 VA € Q x Q, letting successively / <G C£°(M.) such that f(s) = s 



and f(s) = s 2 if \s\ < r y \\X\\g x g, we get that 

(Mt(t),\) g(Xgi = Mi(t) = (Yt(t)-J*Bt(Yt(s))ds,\) g(Xg( 
is a continuous martingale with quadratic variation given by 

< Ml > (t) = [ (A(Yt(s))(\), \) g ^ g da. 



(i 



Proposition 6.7 V/e C b 2 (£ c X <? e ), 

M)(t) = /(Y«(t)) - /*LV(Y«(a))d« 
Jo 

is a continuous martingale, where VY G ^ x CJ^ 

L«/(Y) = ^race (A(Y)£ 2 /(Y)) + (B«(Y), AfOO)^ 
Moreover, the martingale problem associated to the generator L^ is well-posed. 



Proof (of Proposition 6.7) We begin with the following lemma. 
Lemma 6.7 

a : g e x g e — > l+ {g t x ^) , 

where L\ (g^ x g^) is a set of nonnegative operators with finite trace. 
Proof V(Y, A) E (g% x g^) 2 , we have 

(A(Y)(A), A)^ = Re(((L + K)(T)(X), A)^) 

N 



3,1 = 1 

N 2 



E r?« 



TjA/ — T;Aj 



(59) 



with T = Y 1 + iY 2 and A = A 1 + iA 2 . First, V(j, € {1, ... , iV} 2 such that j ^ Z, is nonnegative 
because it is proportional to the power spectral density of Cji at Pi — Pj frequency. Second, the matrix 
r 1 is nonnegative since VX <E C , we have 



1.4 ^ r+oc 7,4 <.+oo 

tXr 1 ^ = y E / E[C„(0)G i (z)]dzA > J l i = - / E[C^(0)^(z)]dz > 

j,l=l •'° 



because it is proportional to the power spectral density of C<^(z) = Y]j C'jj(z)Xj at frequency, and 
with X 3 = Xj/pj, Vj €{!,.. .,N}. Moreover, 



ir aC e(A(Y))=^r j 1 J [(Y J 1 ) 2 + (Y 2 ) 2 ] < -p i ;, ||Y|| 2 4xe£ . 

□ 



, i ie{i,...,iv} 
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Consequently, following the proof of Theorem 4.1.4 in [30 , (59 1 is a martingale. However, and A 
are not bounded functions but this problem can be compensated by the fact that the process Y^(.) 



takes its values in B 



Moreover, from this lemma there exists a linear operator a from x to L2(G^ x G^), which is 
the set of Hilbcrt-Schmidt operators from x to itself, such that A(Y) = a(Y)oa*(Y). According 
to Theorem 3.2.2 and 4.4.1 in [3D], the martingale problem associated to L^ is well-posed because 

\\v(Y)\\<K{N)\\Y\\g tX g t . 



Let us recall that the process Y^(.) is an element of C([0, +oo), (Br y ,g ( xg ( , dB ry ,s 6 x s 5 ))> anc ^ we 
cannot assert that Y^(.) is uniquely determined. In fact, we need to know if its law is supported by 
C([O,+oo),(0 c x0 c ,||.|| ff<xff4 )). Letting 



/(Y) = ||n(£, k 2 ) ® n(£, k 2 )(Y - y)\\l 6xgi , 



where 



n(£, k 2 ) ® n(£, fc 2 ) : 5 £ x ^ 



Y 



n^XY 1 ) 
n(^,fc 2 )(Y 2 ) 



As Y^(.) is a solution on C([0, +oo), (B ryt g ( x g { , ds ry , x o )) °f the martingale associated to L^, we 
get 

E[/(Y« (*))]= Vt>0, 

and therefore II(£, fc 2 )<g>II(£, fc 2 )(Y 5 (.)) = II(£, fc 2 )<g)II(£, k 2 )(Re(y), Im(y)). Consequently, the process 
Y*(.) is strongly continuous since the weak and the strong topologies are the same on M. N . Finally, 
Y^(.) is uniquely characterized as being the unique solution of the martingale prob lem associated to 
and starting from (Re(y),Im(y)), and that concludes the proof of Theorem 



4.1 



6.3 Proof of Theorem I3~2l 

Let Ho = C N x L 2 (0, k 2 ) and y e 7Yo- We begin by showing the tightness of the process (T 4 (.)(y^)) £ , 
which is the unique solution of the martingale problem associated to and starting from j/* = 
II(£, +oo)(y). As the radiating part 11(0, fc 2 )(T^(.)(y^)) of the process T*(.)(y£) is constant equal 
to II(£, k 2 )(y^), to prove the tightness of (T*(.)(y^))j is suffices to show the tightness of the finite- 
dimensional process (II(fc 2 , +oo)(T^(.)(y^)))g. Let E^ be the conditional expectation given the cr- 
algebra cr(T ? (u)(y«), <u<t). Then, Wt > 0, V/i G (0, 1) and Vs e [0, ft,], we have 



||T^+,s)(^)-T«(t)( y «)|| 2 N 

<E«[||Y^(t + S )-Y^(i)|| 2 

N 

< ^Ef[(Y^(i + S )-Yf(t)) 2 



E«[||Y 2 <«(t + ,s)-Y 2 <«(i)|| 



.7=1 
Z=l,2 

iV 



< 



„t+s 2 

E E * [( y t ^/j(Y«(u))^j J + e| [(m| (t + «) - m« (t) 



3=1 
2 = 1,2 



with VY e 5o x 6o, /J 00 = Yj-. Therefore, using that the process T*(.)(j/^) takes its values in 
B Ty ,H v we first get 



L«/j(Y«(«))dt 



< K hf 



and second, 



Ef [(M*,(f + «)-Af*,(f) 



= Ef 



< AfJ, >t+ s - < M), > t 
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with 



< Mji > t - 



L«(/j) 2 (Y*( U )) - 2/j(Y«( U ))tf/j(Y £ (iO)d U . 



Consequently, the process (T ? (.)(?/^))^ is tight on C([0, +00), (Ho, ||.||w„)). Now, to characterize all 
limits of converging subsequences, let us denote by T°(.)(?/) such a limit point. First, for every smooth 
function / on H.q, for every bounded continuous function h, and every sequence < Si < • • • < s n < 
s < i, we have 



E 

Second, 



(t«( Sj )(2/«), 1 < 3 < n) (/(T«(t)(y«)) - /(T%)(y«)) - £ £ 5 /(T«(u)(^))* 



0. 



sup 

TGB r „, ? 



£/(T) - £ 5 /(T) 



<if sup |A^ - A?| 

j£{l,...,JV} 



I 3 3 1 



I K | K i I 



Consequently, T°(.)(?/) is a solution of the martingale problem associated to £ and starting from y. 
However, following the proof of the uniqueness in Theorem |4.1[ this martingale problem is well-posed 
and therefore T^(. )(?/£) converges in distribution to the unique solution of the martingale problem 
associated to £ and starting from y. 

6.4 Proof of Theorem 15.21 

The proof of this theorem follows ideas developed in j^Zl Chapter 11]. In order to prove this theorem 
we use a probabilistic representation of z) by using the Feynman-Kac formula. To this end, we 

introduce the jump Markov process (X t Ar ) t>Q with state space ( — (N — 1) /TV, . . . , 0, . . . , (N — 1) /N} 
and generator given by 



£ 



N , 



= r 



u+i 



1 - 1 

N 



1+2 1+1 



for I € {!,..., N -2}, 



£ 



N , 



— r c 

— 1 1, 



|2|+2|Z| + 1 

farle{-(JV-2),...,-l}, 
and 



-N , 



N 



I - 1 
N 



N-1N 



\l\\l\ + l 



0(0) 



±(N- 2) 
TV 



A 



1 + 1 

JV 



Z + l 

AT 

0(0) 

±QV-1) 
AT 



Using the Feynman-Kac formula, we get for (j, I) £ {1, . . . , A(o>)} 2 



- A N Jo 1 (|X«| = ^=li) dl ' _A «-l ■'o' 1 (| 



\ dv 



1 (l^l+^ = ^) 



Let / be a bounded continuous function on [0, 1], we consider T l (us, L) as a family of bounded measures 
on [0, 1] by setting 



7j(w,L) = Ei 



A N Jo 1 ^ x N ]= N__l^ dv - k N-l Jo 1 



(^^)*/(|Xf| + i 



In the first part of the proof, we consider the case v £ [0,1) and in a second part we shall treat 
the case v = 1. 
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Let u £ [0,1) such that l(N)/N — > u. We begin by introducing some notations. Throughout 
the proof we denote by Tj) N the Ith passage in j/N, for j £ {—(N — 1), . . . , TV — 1}. To avoid the 
unboundness in C N of the reflecting barriers C N <f)(±(N — Vj/N), we introduce the stopping time 

r a - r (1) A r (1) 

N — '(N-[N a ])/N n -(N-[N°>])/N 

with a £ (0, 1). Let Zf' T = X& T c, Vi > 0, be the stopped process and d(N) = (Z(iV) - l)/iV. We 
denote by ^d(N) the law of (X^) t >o starting from d(N) and by P^n the law of (X^' T ) t >o starting 
from d(N). Let 



dv \ dv 

where a,oo{-) £ C 1 (]R) is an extension over R of aoo(-), which is defined on [—1, 1], and such that the 
martingale problem associated to and starting from u is well posed. We denote by P„ this unique 
solution. Let ip £ C^(R), 

M v (t) = (fi(x(t)) - <p(x(0)) - f Ca^{x{s))d Sl 

Jo 

and r r = inf(u > 0, \x(t) \ > r) for r £ (0, 1). 
Lemma 6.8 Vip £ Cg°(R), Va £ (2/3, 1), 

lira sup \£ N <f(v) - £a x (p(v)\ = 0, 

ut [ JV > N J u [ ]V > N J 

where £ N ip(v) is defined as follows. Vj £ {1, . . . , N — 2}, 



■ r j-+ij+2 M ^xT- ) -V 



at y m at 



/or « e [j/AT, (j + 1)/N), and 

+ ^+l i+ 2 (*>( 

ifve (-(j + l)/N,-j/N]. 



-j - 1 \ f-j 



N J \ N 



Proof (of Lemma 6.8) We shall restrict the proof of this lemma to the proof of 

lim sup \C <fi(v) — C T i <f{v)\ = 0, 



since the other case is completely similar by symmetry. We start the proof of this technical lemma 
by proving Lemma |2.1| page [5] 



Proof (of Lemma 2.1) Let h e (v) — . v == and g(v) = arctan(w). We recall that Vj £ 

' ' yj (nikdB) —v 2 

{l, . . . , AT}, tan(CTj) = —h e (<jj) First, 

\a j+1 - o-j -7r| < |g(tan (a j+1 - {j + 1)tt)) - 5 (tan (cr,- — J7r) ) | 

< AT|tan (aj+i) — tan (<Jj)\ 

< K\h e (a J+ i) - h e (<Tj)\ 

< K sup h' e (v) 

»6[(Tj,ITj + J 
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where h'Jv) — {nikde) — which is a positive and increasing function. Moreover, 

eV ' ({ ni kd8) 2 -v 2 )2 



and then 



sup 

je{i,...,N-[N<*]} 



cr N _ [Na] <(N-[N a ])7T 



Second, in the same way we have 

<Tj +2 -2<Tj + l + <Jj 

= ,g(tan (a j+2 ~ (j + 2)tt)) - 2g(tan fa +1 - (j + 1)tt)) + 5 (tan fa - jir)) 
= -(g{he{o-j+2)) - 2g(h e fa+i)) + g(h e (<jj))) 



and 



g(Kfa +2 )) - 2g(h e fa +1 )) + g{h e {a 3 )) = [g'{h e fa +l )) - g'(h e fa))] . [h e (a j+2 ) - h e (a j+1 )] 

+ 9'{h e ((Tj)) [h e fa +2 ) - 2h e (cr j+1 ) + h e fa)] 

(h e fa +2 )-t)g"(t) dt 



M°j+i) 

h e (<Tj+l) 
h e {tTj) 



(h e fa +1 )-t)g"(t) dt. 



Moreover, 



\g'(hefa+i))-g'(h e fa))\]h e fa + 2)-hefa+i)\<KN 1 - 3a . 
hefa+z) - 2h e (a j+1 ) + h e fa) - 



with 



a-j+2 



h' e (t)-K(t-n)dt 



a ' +a ti B (t)-h' e (t-n)dt 

o-j + l 

+ / h' e (t - 7r) dt - / h' e (t)dt, 
<Kh'l{cx N _ [Na] ) = 0{N 1 *-i a ), 



because h"{v) 



(( ni kd0) 2 -v 2 )5/ 2 



rcrj+2 raj + i+ir 

j h' e {t - it) dt - / h' e (t-7r)dt 



Finally, 



< h' e (a N _ [Na] ) (\<T j+2 - (T j+ i - tt| + \a j+2 - <Tj+i - 7r] 



h e ((T J + i) 



(fte(o-j+2) - t) g"{t) dt - / (hefa+l) - t) ff"(t) dt 



< 



K ((hefa+2) - h e (cTj + i)) 2 + (h e ((T j+1 ) - h e ((Jj)) 2 ) 



l-3a 



and 



< KN 



sup \<t j+2 - 2a j+1 + <Tj\ = 0(N 1 ~ 3a ). 

j£{l,...,N-[N°>}-2} 



This completes the proof of Lemma 2.1 since we can take a > 1/3 and we have N2~2 a < TV 1 - 3 ". □ 
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From this lemma, we immediately get 



sup 

je{i,...,iV-[Ar<*]-i} 



Before showing that 



\S((T J + l - <Jj,(Jj + l - (Jj) - S(tT,Tt)\ 

< K sup 

je{i,...,iv-[iv<»]-i} 



\a j+ i -dj — tt| =o(N 1 i-i a ^j 



sup 

where Aj is defined by pi, we prove that 



A i d 



OiN 01 - 1 ), 



SUp \(Tj — J7T | = O 
j£{l,...,[iV«]} 



1 



N l-a 



In fact, Vj G {l,...,iV Q } 



\<jj — J 7T | — |g(tan {(jj — J7r)) — arctan(tan(0)) 
< if| tan((Tj) | 



< 



Kh e (a[ Na] ). 



Moreover, <J[n<*} < N a ir and then 



fc (^])=°(]^) 



Consequently, 



because 



sup 

je{i,...,N-[N°]} 



sup 

je{i,...,N-[N«]} 



A'- 2 



< K sup 

je{i,...,N-[N"]} 



sm 2 (<jj) sin(2(jj) 







2(7, 



< 



A' 



sin 2 (aj) sin(2(7j) 







2(7, 



< 



(mkd8) 2 - <J 2 N _ [Na] 



< 



sup 

j6{[JV<*] + l,...,iV-[JV<*]} 

K 1 



sin(2(Tj) 



2(7," 



sup 

36{1,...,[JV«]} 



sin(2(7j 



2(7 7 



N l/2+a/2 2(7 



sup 



< if 



[jv<*]+i je{i,...,[iv»]} 
1 1 



sin(2<Tj) — sin(2j7r) 



2(7, 



iv^ + ^ + ^, e{1 ; u ^ ]} '^- j7r ' 



Now, let us introduce 



D 



N 



\A 2 jA 2 j +1 S(cij + i - <Jj,a J+1 - Oj) 



2n lV 1 nf^d 2 A/ 1 .,;;,-:,/•■ 



a + 08,- - /3 J+1 r 



Then, for j G {l, . . . , - 2} and v e *±i] 

2 / TVtt \ 2 



jv fnikd6\ 2 [ Ntt\\( f[Nv\ + l\ (WvU\ N 



D 



N 
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Consequently, from the following decomposition 



JV 



TV 



JV 



f +1 + (K 



(T)-M^)^(^) 



[JVv] 



TV 2 

■at - , , - , 

N I V JV 
n l d202 f[Nv] + l 



[Nv] - 1 
JV 



■ [Nv] - 1 > 



n 2 d 2 e< 



C aoo ip(v) 



[Nv] + 1 \ / [JVj)l 
^l 1 tr l-W 



N 

N r>N\ 



nfd 



B N _ ■_^_ aoc(v) 



n\d 2 d 2 d 



l^ aoo (v)[N 2 ^ 
(v)[n(ip( 



N ( B j+i - B j) ~ n2 dy 
[Nv] \ ( [Nv] - 1 



a-oo(v) 



n \d 2 e 2 d 

7r 2 dv 



N 

[Nv] + 1 
JV 



N 



fWv]\\ , 



sup 

' 1 N-[N a ] - 
N ' N 



N 



and because it is easy to show that 

' [Nv] + 1 
_ N 

[Nv] + 1 

N ' 



sup 

j_ JV-[JV° 
JV ' N 



= o 



N 



[Nv] - 1 

N~ ' 



<p"(v) 



O 



it suffices to show the two following points 



lim 

JV 



sup 



sup 



je{i,...,iV-[A™]-i} i+i] 



D 



N 



n\d 2 6 2 



■aoo(v) 



lim 

N 



sup 



sup 



je{i,...,N-[N°]-i} ve \±m.] 



N[Bf-Bf +l 



njdW d 

7T 2 dv 



aoo{v) 



= 0. 



We decompose the proof of these two points into two sublemmas. 
Lemma 6.9 



lim 

N 



sup 



sup 



je{i,...,;v-[iv«]-iwi i+il 



n n 2 d 2 e 2 , . 



= 



Proof (of Lemma [679]) Vj € {l, . . . , JV - [JV Q ] - l} and £ 
inequalities, 



J_ 2+1 
AT ' AT 



we have the following 



< 



2 k 2 d 2 



< 



n 2 k 2 d 2 



cr 2 ,, \ 1 _ + 1 



Moreover, for I £ {j,j+ 1} 



1 - (9vy 



1 - 



< 



n 2 k 2 d 2 



riikd 



N 



26 



{\-e 2 y 



< 



K 
JV' 



Consequently, 
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Next, 



n\k 2 d 2 



nfk 2 d 2 J d 



n\kd 



2 k 2 d 2 



< K 



1 



'3+1 



K 



< 



K 



n\k 2 d 2 

(<7j+i - (Tj - tt) 

1 



i\k 2 d 2 ) d 



n\kd 



n'ik 2 d 2 



n\kd 



and then 



^(1-02)3/2^ + 1-^ 
1 



n 2 k 2 d 2 
2 + K Ci'+l — (Tj — TT 



sup 

je{i,...,N-[N<*]-i} 



1 



a 2 + (ft - j+1 y 



,2 I 

1 + d 2 



i 2 k 2 d 2 



Moreover Vj € {l, . . . ,iV — [7V Q ] - 1} and Vu e 



j j'+i 

AT' TV 



™ 2 fc 2 d 2 

we have 



n\kd 



Ov 



1 - 



< if 



(T, 



n\kd 



6v 



< 



K 

N' 



n 2 k 2 d 2 



and finally 



sup sup 

jG{l,...,N-[iV°]-l} [i i+l 



^j-n/ a -t- rf2 1 _( 0a . ) 



= O 



(jvs-S«). 



This concludes the proof of Lemma 6.9 □ 
Lemma 6.10 



lim 

N 



sup 



sup 



NiBy-B?^ nld2 ° 2d 



J j+i 



TT 2 dv 



aoc(v) 



Proof (of Lemma 6.10) We separate the proof of this lemma into two step. First, for each j 
{!,..., N - [N a ] -2pet 



Cf = N 



\ 



k 2 d 2 



n'fk 2 d 2 



3 , / I _ J + l 

n 2 k 2 d 2 y 1 n 2 fc 2 d 2 / 



We can write Vw S 

2<9 2 v 



C 



jv 



3 (1 - (0v) 2 ) 2 



,N 



N 



1 



% 2 k 2 d 2 



"3 + 1 



' T j I - 
ni kd 



dw — 



Jfa (l-w 2 )i~~ (l-(9v) 2 ) 2 



kd 



-dw - N 



\m kd 



Ov 



N 



(l- w 2)f \ ni kd n x kd) (l _ (9v) 2 )i 

Ov / 1 1 



N 



n\kd n\kd) (1 

°J+2 _ Uj 

nikd ri\kd 



\ Ov / 1 

/ n - r^^l v / 



1 ri 2 /c 2 <i 2 



26* 



<9r 



(1 - {6v) 2 Y 
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We can check that the function v i— » — is bounded on [0, 1] and 



(l-(8v) 2 ) 2 

N 



< 



Moreover, v 



n\kd 

is bounded on [0,1] and 
1 1 



Vj+2 - (?j - 2?r 



+ 29 



Nit 



nidkO 



- 1 



(i-(<M 2 ) 



Finally, < 



N 



< 



y/1 - (9v) 2 
and 



< 



k e 



N (1 — 6)2)3/2 - 



0v 



" 1 ''''' 



I \nikd mkdJ (1 - (0u) 2 )§ 



< AT 



2# 2 



(l-02)§ 

( mfcrf) ( nikd) 



< 



N 



Consequently, 



sup sup 

ie{i,...,*-[JV«]-i}„ e [^] 



' (1 - (^) 2 ) 2 



= 0(7Vs-i a ). 



Second, for j e {l,...,N- [N a ] - 1} and v € [jf, ^] , we have 
1 1 



N 



7T 2e z -» 
d 2 (l-(0i;) 2 ) 2 



a 2 + (ft+i - ft+2) 2 a 2 + (ft -ft+i) 2 / / 2 7r 2 (<fa) 2 ^ 

^ M d 2 i-(e«) 2 y 



< 



N 



(a 2 



1 



+ 0ft +1 -/3 j+2 ) 2 a 2 + (ft - 



- ft +2 ) - (ft - ft + l)) , 2 ^ 

(a 2 + - ft+i) 2 ) 



+ 



M(ft +1 - ft +2 ) - (ft - ft +1 )) . - 2 ^-^„ + 5 (1 " ;;)2)2 



(a 2 + (ft-ft +1 ) 2 ) ( a 2 + f^ 



For the first term on the right of the previous inequality, we have 
1 1 



N 



(a 2 



+ - p j+2 y a 2 + {Pi j+1 y 



-2(/J,- - /3 j+1 ) 



(a 2 + (^-/3, +1 ) 2 )' 
<KN({P 3+1 -P 0+2 )-{(3 -(i J+l ))\ 
and we shall see just below that 

sup \(3 j+2 -2(3 j+1 -(3j\=o(^). 



je{l,...,N-[N°]-2} 

Now, for the second term we have previously get 



sup sup 

,e{i,...,AMA-MU e [^,+i; 



ft ~ 03+i - 



7T 8v 



= 0(ATs-i a ). 
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Then, to finish the proof of this lemma it suffices to show that 



sup sup 



(1 - (9v) 2 ) 2 

To show this relation we shall use the following decompositions. For I £ {j,j + 1} 

1 



= 0(N 2 ~ 3a ). 



'1 



i 2 k 2 d? 



1 



'i+i 



n\k 2 d? riikd 



(a i+1 - ai) 



n\kd 



1 - 



"i+i , 



\n\kd 



w 



n\k 2 d 2 
1 

(1 -w 2 ) 



-dw, 



and 



Nnxk(\\l- 
N 

: 1l 



i\k 2 d 2 



-2W1- 



"j+i 

?l 2 /fc 2 d 2 



(Vj+i -o-j) 



riikd 



V n 2 /c 2 d 2 / 



(1 - (6»w) 2 ) i 



°J + 1 

nikd 



Wd 2 



+ Nruk 



r j+i 



\nikd J (1 



n'fk 2 d 2 

7 3 +2 

"»i*<l / (Tj-\-2 



(1 - (6>?;) 2 ) ! 
1 



i±i Vriifcd / (1 - w 2) 

i fed * ' 



First, using Lemma [2. 1| we have 



nikd 



w 



(i - w 2 )-; 



"J+2 



(1 — If 2 ) 2 



Second, we have 

"(Vj + l - CTj 



7i i k 



1 - 



1 _ "j+i 
1 n 2 k 2 



(1 - (^) 2 ) 



2\2 



N 



n\kd 



mkd 



Wd? 



1 - 



N 



(°J+2 - 2CTj-+1 + o-j-) 



7i( N 



nikd 



CT j + i 
n\kd 



2 k 2 d 2 



(1 - (6»w) 2 )f 



where, according to Lemma 2.1 the first and the third term are 0(N 2 and the second term 

0(N 2 - 3a ). That concludes the proof of Lemma |6T0] for a € (2/3, 1).D 

Consequently, thanks to Lemma |6.9| and Lemma |6.10| we get 



sup 

i i 



\C N <p(v) - £ a ^(v)\ = (N( 2 -^( a -i)^ 



this concludes the proof of Lemma [6.8| D 
Lemma 6.11 is tight on X>([0, +oo), 
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Proof (of Lemma 6.11) Let M. t = a(x(u), < u < t). According to Theorem 3 in [TBI Chapter 
3], we have to show the two following points. First, 



lim lim 



sup|a;(t)| > K = 



The first point is satisfied since we have VAT, P^'^ (sup t>0 |:z;(i)| < l) = 1. Second, for each N, 
he (0,1), se [0,h] andi>0, 

E r ^((x(t + s) - x{t)) 2 \M t ) <Kh. 
Concerning the second point, letting (p E C£°(R) such that tp(s) = s if \s\ < 1, we have 

E P ™<«> ((x(t + s)- x{t)) 2 \M t ) < 2E P ~<«> ((M*(t + s) - \t)) 2 \M t ) 



2EV) 



t+s . - 

C N ip{x{w))dw 



M t 



with 



N 

d(N) 



which is a (A^t)*>o-martingale under P^(N) anc ^ we know that 



Pj^sup^/.j ; 



N - [N a ] 
N 



1. 



Moreover, by Lemma |6.8| 



sup sup 

" e N • N \ [N ' N 



\C ip(v)\< +00 



and the fact that C y(0) = 0, we get 



E r d(«) 



t+s 



£ N tp(x(w))dw 



M t < Ch 2 



We recall that 



< > t = / (C N V 2 - 2<pC N V ) (x(s))d. 



Then, using the martingale property of (M^ (i))t>o, we have 

E P ™<«> ((M*(t + s)- M^(t)) 2 \M t ) = E p V> {{M${{t + s) A r£) - Af *(t A T%)) 2 \M t ) 

= E p "<«) (M*((i + s) A r£) 2 - M"(t A r^) 2 \M t ) 



E r -<«) (< M^ v > (t+s) ^ N - < M» > tAT « \M t 



= E rd <«) 
< Ch. 



tAr° 



In fact, by Lemma |6.8| we have 



sup sup 

N r r jv-[jv°] i 1 , i f i Jv-[«°] 1 



<+oo, 



sup 



sup 



|£>>)| <+oo, 



in addition to C N ip(0) = and swp N C N (p 2 (0) = & < +oo. □ 
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Lemma 6.12 Let Q u be a limit point of the relatively compact sequence (^^^^j ■ Then, \f(p € 
C^°(1R) and W G (0,1), (M v (t A r r )) t > is a (M)t-martingale under Q u . 

Proof (of Lemma 6.12) Let (v^^,^j be a converging subsequence. Throughout this proof we 

will take N for N' to simplify the notations. Let < t\ < and $ be a bounded continuous 
.M tl -measurable function. We have 



Furthermore, Vi > 



E r d(N) 



E 



£V(»))l(x(.)ejj}) 



ds<I> 



tAi> 



£ ^(a;(s))l (a: ( s ) = o)ds$ , 



with 7^ = [-(N - [iV Q ])/AT, -1/iV] U [1/N, (N - [N a ])/N]. Using Lemma[6. 

ctAr r 



lim 



P JV,T 

E *M 



Consequently, we have to prove the two following points: 



= 0. 



lim 7V E f V«> (MJt A T r )$) = E^" (M„(* A T r )$). 



. li mj vE P "<«) (/ * Ar " l (x(s)=0) ds) = 0. 

We prove the first point as follows. The problem is to apply the mapping theorem to the functional 
M v (t/\T r ) and to do this we must have Q«(-Dm (tAr r )) — 0, where f m (tAr r ) is the set of discontinuities 
of M v (t A r r ) for the Skorokhod topology. While M v (t) is continuous for this topology, it is not 
necessarily true for t t . However, we can follow the proof of Lemma 11.1.3 in [27 and then use a 
family of stopping times for which we can apply the mapping theorem. 

We know that the size of the jumps of (X^) t is constant equal to 1/N, therefore we have 
Q„(C([0,+oo),R)) = 1 (see Theorem 13.4 in [3 for instance). Then 

Qu(D Mv{tATr) ) = Q„(£> M „(tAr r ) nC([0,+oo),R)). 

We recall that the Skorokhod topology on C([0, +oo),R) coincides with the usual topology defined on 
this space. Therefore, £*M v (tAT r ) H C([0, +oo), R) is the set of discontinuities of M v (t A r r ) under the 
topology of C([0, +00), R), and r r restrict to C([0,+oo),R) is lower semi-continuous. Consequently, 
according to the proof of lemmas 11.1.2 in (27], there exists a sequence (r n ) n such that r„/r and 

Q«((r rn <+oo) r\D Trn nC([0,+oo),R)) =Q„((r r „ < +oo) n D Tr )) = 0. 

Then, Q. U (D M ttAr r )) = an d we can a Pply the mapping theorem to M v (t A T rn ), i.e 

limE P "(«) (M v (t A r r J$) = E Q " (M v (t A r r J$) . 

Finally, we obtain 

E Q « (M v (t 2 A 7- r J$) = E Q " (Af v (fi A r r J$) , 

and 

limE Q " (MJt A r r „)$) = E Q " {MJt A r r )$) 

n 

because T Tn f r r . Consequently, 



E Q " {M v {t 2 A r r )$) = E Q " (M v (ti A r r )$) 
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For the second point, we have 

ft/\T r 



E 



™ Of 



nt/\T r 

Jo 



< En 



since the stopped process spends less time in than the original process and the last inequality is 
given by the Markov property. We denote by JV t the number of returns in during the time interval 
[0, t] and by (5^)j>o the renewal process associated with the return times in 0, (cto'^)»>1) of the process 
(X t N ) t , with Y = c^ 0) = 0. Moreover, for a' e (0, 1) 



E 



f l W=0) da <tP (jV t °> [N 1+a ']\ +E £ r 

•'0 J ,_n "^fn 3 



[N 1 + a ] ,_« + !) 



L(X s «=0) 



< 



^ Eo M + |A " + "' l + 1 



[jv 1 +«']" ul " tJ 1 r§! 

since ( / ° j} l/ X N =o)ds) is an i.i.d sequence with mean l/rs^. We recall that iV t ° + 1 is a stopping 

\ CT o 3 '3 

time for (lj)j>i. Then, 



Furthermore, 



En 







E n 

E 
E 



o 



= En 



N° + l 



(E [iV°] + 1) E 



S 1 ) 



J N ? +1 ) ft I 1 

mf(s>T t N ,X? +s = 0)l( x „^ 



+ E 



inf( S >0,^ s = 0)l (X(Vo) 



where T^ = inf (s > 0, X^ +s ^ 0) . Then, using the Markov property we get 



En 



and 



En 



inf( S >0,X^ s = 0)l (X(Vo) 



xf=o)J = ( t - 


-E 


[< 


)Po(xf 




\ 27V - ] 

= 0)+ T e 
1 12 


N-l 








= E E " 


inf(*>0,X& a = 


j=-(AT-l) 








N-l 








■ E 


FE j 


T (1)_ 


)Po(Xf 


J=-(AT-1) 









JV j 



EE^ p o(i^i = i)+*Po(xf^o) 



r c 



+ iP (Xf ^0) 



where if is a constant independent of AT. Consequently, 



and 



□ 



E [7V t °] < A" 

i 



1 1: 



2AT- 1 



E 



2JV-1 rj 2 



1 



jya'A(l-a') 
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From Lemma 



6.12 



we have Vr G (0, 1), Q„ = F u on A4 Tr . From this relation and the fact that 
Q u (C([0,+oo),K))^P u (C([0,+oo),K)) = 1, Q u = ¥ u on M T [ since r r / n as r / 1. 

Let / G C°°([0,1]) with compact support included in [0,1) and let (^^ffi^j be a converging 



subsequence as in the previous proof. We have 



7? N '\u>,t)=E d(N , ) 



f \x 



N' 



N' ) ( t<T N>) 



(60) 



with 



r(N) = E d{N) 



-A«/o l( |x ™ | = iv-i) 



x / ( | 



N 



N = T (N-1)/N A T -(N-1)/N and A G Usin § Lemma 



where t\ 

the first term on the right in (|60|) 



\T<*<t<T%+\) + !(t>T^+A) 

and Lemma 



6.11 



6.12 



we can study 



E d < N ') 



J.A 



since (t < t%) G M r; and * d{N) - * d(JV) 



on A^ r ° . Moreover, 



/( 1^)1 + ^7 



= E<"0 [/(|x(t)|)] + (l) 



-/ 



AT' - [iV'«] + 1 

W' 



Consequently, 



However, 



limE 

iV' 



d(iV') 



= E^[/(|z(i)|)]. 



[/(|i(t)|)l(n<t)] =0. 



In fact, let r s = inf(t > 0,Vi> > t,x(v) = x(t)) be the first time for which the process becomes 
constant. From the Portmanteau theorem 



1 = lim 

N , d(N 



i?) ((r. < n)) < Q„ ((r s < ri)) 



where A denote the closure under the Skorokhod topology of a subset A of 2?([0, +oo),K). Moreover, 
we have 

(t s < ri) n (n < t) n (x(o) = u) nC([o,+oo),m) 

= (r s < n) n (n < t) n (x(o) = u) nC([o,+oo),R). 



Then, 



K(\x(t)\ e supp{f),T X <t)< Q u (\x(t)\ e su PP (f),T s <n<t) 

<Q u (\x(t)\ G S? w(/),K*)I = 1) =0, 



and 



limE, 



AT' J x (*<^) 



= E^ [/(K*)|)l (t < Tl) ] = E f " [/(|ar(t)|)l (t<7l) ] 
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Finally, by the following lemma we get 



lkaT« N '\u,,t)=E W - [/(>(t)|)l (t<Tl) ] 



We can remark that this limit does not depend on the subsequence (N 1 ). The following lemma 
represents the loss of energy from the propagating modes produced by the coupling between the 
propagating and the radiating modes. Moreover, this lemma implies the absorbing condition at the 



boundary 1 in Theorem |5.2| which implies the dissipation behavior in Theorem 5.3 

Lemma 6.13 lim N * r(N') = 0. 

Proof 



|r(JV")| < H/lloo E d(JV0 



~ A «< /o 1 



> ds 



(t>r° +A) 



- d(N') 



\X t N '\ + supp(f),r%, <t<T° N , + X^ 



First, let a' S (3/4, 1) and 7V t the number of passages in (N — l)/N during the time interval [0,t]. 



E. 



d(N') 



~ A N' So 1 (\ X M' l= N'-l\ ds 



(t>T°,+A) 



< E 



d(N') 



~ A «< So 1 



ds 



X |1 (*>^_ 1) ^+A) + 1 ( tL; 



(1) 

(N' -1)/N 



,+A) 



We shall work only with E^jv 
for the other term. 



L (t>^_ 1)/N ,+A) 



but the same proof works 



E 



d(N') 



~ A N' So 1 (l X N'l= N '- 1 \ dS -. 
P, \ 1 s N' } 1 / , •> 



< E 



d(iV') 



(N'—l)/N' 
~ A N' Jo 1 (' X N' = «'-l S | rf;i -, 

e 1 3 ^ j J-CATf >[AT-']+l) 



+ P d(A r,) (iVf < [JV a '], 



1 r (JV'-l)/iV' - A 



On (iVf > [iV' Q ] + 1), we have 



JV' ^ + 1 > 



[AT' a 'l AC f '(JV'-l)/JV' 



,0+1) 



< n e 



(]Y'-1)/W V s N / 



We denote by the time of the first return in (A^ — 1)/^, then 



_A N' So 1 ( X N' = N'-I\ ds 



IN'" ] 

< n e<hn>) 



(iyf'>[JV"»'] + l) 

,c j. (N'-l)/N' . 
A N'J U) (xN' = N'-l \ 

(N'-l\/N> V s N ' 



< E iV /_ 1 



A c r {N'-l)/N' 
^JV' JO |X»'= N M 

3 V s N' } 
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since f , < ™ _1)/JV 1/ YN _ N-i \ds I is an i.i.d sequence. Moreover, we can check that A c „ > CN 3 ' 2 and 



then 



E iv'-i 



.c (• (JV'-1)/N' , 
~ A N' Jo 1 ^ X W' = AT'- 1 ~) fflS 



< Ejv'-l 



(1) 

13/2 r"(N'-l)/N' 



-CN' 3 ^J 



In fact, a computation gives 

afc 4 A 2 mikd 



A c — 



N 



Vy^V 2 ~ (nikdO) 7 



16n(3 N J nikdg ^( ni kdy-ri 2 (i + (/fo - ^(m/firf) 2 - r? 2 ) 2 ) 

^(jy - <tn,V- &n) 



(j] 2 — (nikdO) 2 ) sin (77) + rj 2 cos 2 (r/) 



drj. 



However, we recall that the support of S lies in the square [— , ] x [— ^f, 4f] , then we can restrict 
the integration over \n\kd6, n\kd9 + Moreover, 



[P N - \J (nifc) 2 - r) 2 /d 2 ) = (tj - cr) 2 ^ - ^ 2 , for some y e 



(Tat ?7 



nifcci' n\kd 



'3iA 2 V 2 /( ni kd) 2 
2 J l- V 2 /( ni kd) 2 ~ ' 



where K stands for a constant independent of N, because 9 < 1 and fc 3> 1. Therefore, 



A° N > K 



7 

J n\kdO 



v 



vV - {mkdeydr) > K"N 3 / 2 , 



where we assume that the function S has a positive minimum, and then K" > 0. 
Now, let us remark that Vw e [0, h^iV 1 / 4 )] 



e - " < 1 



ATi/4 



and 



Then, we get 



A" 



(1) 

f a {N' -1)/N' 

Jo V'=^: 



1 N'-1N' 



Ejv'-i 



><-// M l3/2 r (N>-l)/N> . 
~ K " JO 1 ( V N>- N'-l \ as 

e \ x ° -Tv^J 



if" 

< 1 TTT [ 1 



and 



E d(iV') 



- A JV' Jo' 1 /'^]v'_ iV'-l \ ds 

e l x = -^v^J 1 



(ATf >[N'"'] + 1) 



< e 



iV' 3 / 4 V In(JV'V4) 



Moreover, 



d(JV) 



(< < [A^],*- r$_ 1)/JV > A) < P d(JV) 



_(["«]) _ (i) >A 



<r 1 TB> f ([W'l) (!) 

- ^^(JV) I T (Ar_i)/Ar - T (AT-l)/W 



7V C 



< — —^(N-l)/N 



r (1) 

'(N-l)/N 



< 



K 

7V 1_a ' ' 
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Consequently, 



limE 

N> 



d(N') 



, ds 



A N> So 1 (\ X N / I- N'-l \ 



= 0. 



Second, let c/ G (0, 1) such that supp(f) C [0, c/ — 1/iV'] and x £ [0, c/), then 

P d{N) (\X t N \ + j^£ SU PP (f),T% <t<T° N + X S j< ¥ d{N) (X t N £ [-Cf,Cf],T% <t<r° N + X) 



< 



- d(N) 
"d(JV) 



X? €[-c f ,c f ],X? 



N 



N — [N a ] _ Q 



N 



') 'N 



<t<T% + \ 



X? £ [-c f ,c f ],X» = - - } Na] ,r% <t<T% + \ 



N 



We shall treat only the case where X^L = (N — [N a ])/N, but the following proof works also in the 
other case. Let c/ £ (c/, 1), p £ (0, 1) such that [5/ — p, c/ + p] C (c/, 1) and A' £ (0, 1). Using the 
strong Markov property we have 



d(N) X t £ [-C/,C/],X r o = 



N 



,T% <t<T^ + X 



< 



^ (t$-i),n > A ') + P %/I ( t ^/±p < A + A % 
where Tc f ± p — inf(i > 0, \x(t) — Cf \ > p). First, a computation gives 



1 

A 7 



r (1) 

'(N-l)/N 



_ 1_ y-? A^+l + l X 



l=N-[N a ] i + li+2 



Second, the sequence (r(JV'))jy/ is bounded. Let (r(iV"))jv» be a con vergin g subsequence. We recall 



on -M r =, where c(iV) = [iVc^]/JV, and by Lemma 



6.11 



the sequence (p^jv")) 

is tight. Let (P^iV'") ) be a converging subsequence to Qg / . Moreover, Tc f ± p < and therefore 



that IP / \r\ — IP / nr\ 

c[N) c(N) 



( T c/±p < A + A') £ .M. T ". Consequently, by the Portmanteau theorem 



lim 

N"> 



j AT'" 



P') K±p < A + A') = lim P^V (r £/±p < A + A') 



< lim 



< 



i^/'/T) ((t£/±p < A + A')) 
((r £/±p <A + A')). 



We recall that Qg, (C([0, +oo),R)) = 1 and we can show that 



{t~c s ± p < A + A') nC([0,+oo),R) = (t £/±p < A + A') n C([0, +oo),R). 



Then, 



and 



%%V$m>) fe/±? ^ A ') ^ % fc/±/> < A + A') 



Hmr-(iV"') < Q Sf (r~ Cf±p < A + A') . 
Finally, limjv'» r(N"') = and the limit of all subsequences (r(iV"))jv" of (r(N')) N * is 0. □ 

To finish, is a bounded sequence. Let (t^ N \uj,t)J be a converging subse- 

quence. By the previous work, there exists an another subsequence such that 

fimr; (JV '' ) ( W ,<) = E f "[/(| 2 :(t)|)l (t<ri) ], 
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where the limit does not depend on the particular subsequence, then all subsequence limits of 
uJ,t)j are equal to E p " [/(|ar(i)|)l( t<T1 )] . Consequently, 



limT; (JV) (^)=E p " [f(\x(t)\)l {t<Tl) ] 



N f 

Now, we have to show that this equality holds even for a sequence (1(N))n such that l(N)/N 
f 

(i) 



u ~ 1, i.c limjv T 1 S N \lo 1 t) = 0. To do this, we write for A G (0,t) 



r; (A Vt) < ||/||oo(p dW (t < r ( ( ^ 1)/JV) + a) 



•E 



d(N) 



-A N J* 1/ JV _ JV _ 1 N(iu 



( t >^i)/« + A J 



We have already shown in Lemma 6.13 that the second term on the right in the previous inequality 
goes to 0. The proof did not depend of the sequence (d(N)) n. Moreover, we have 



A < 



1 



-E 



t- A 



d(N) 



r (1) 



and 



E 



d(N) 



r (1) 

r (N-l)/N) 



N-2 



Consequently, we have Vu G [0, 1] and V(1(N))n such that l(N)/N — > u, 

limT; (Ar W) = E'* [/(|^)|)l (t<Tl) ] , 



(61) 



where the limit satisfies the required conditions. Finally, from the decomposition used in the proof of 
5.3 we have V</? G L 2 (0, 1) and (p a smooth function with compact support 



Theorem 



\\T»(L, .) - %(L, .)|| L2(0 ,i) <2\\<p- 011^(0,1) + ||T/(L, .) - 7$>(L, .)||l»(o,i). 



Using the density of the smooth functions with compact support in L 2 (0, 1) for ||.||l 2 (o,i) an d the 
dominated convergence theorem we get the first point of Theorem |5.2| The second point is a direct 



consequence of the probabilistic representation (61 1 and the density for the sup norm over [0, 1] in 
{f G C°([0, 1]), ip(l) = 0} of the smooth functions with compact support included in [0, 1). 



6.5 Proof of Theorem 15.41 

As in the proof of Theorem . 



5.2 



we use a probabilistic representation of T^'\z) by using the Feynman- 
Kac formula. However, we introduce the jump Markov process which is a symmetric version with 
respect to reflecting barrier [N — 1)/N of that used in the proof of Theorem |5.2| 



Let (X^) t>Q be a jump Markov process with state space { - (N- 1)/N, . . . , (JV- 1)/N, . . . , 3(N 



1)/AT} and generator given by 



C 



= r 



|/|+2|Z| + 1 



I - 1 

N 



forZG{-(iV-2),...,-l}, 



for le{l,...,N-2} 



= r 



u+i 



i - i 

N 



\i\\i\+i 



1+2 l + l 



= r 



|i-2(JV-l)|+2 |i-2(W-l)|+l 
1 |J-2(JV-1)| |i—3(JV— 



I - 1 

N 
l + l 
N 



l + l 

N 



l+l 

N 
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for I e {N, . . . , 2JV — 3}, 



. ( i 



r ?-2(W-l) l + l-2(N-l) IV I I - V I ]y 



for Z G {27V- 1,..., 37V- 2}, 



Al-^Ur ^ N - 1 



N I -N-iNyry N I T \ N 



3JV-3 \ / / 3JV-4 \ J 3N-3 

l ' ;| t n^v- -n^v- 



^(o) = ^(,(^)-,(o)) + ^(,(^)-0(o) 

V - ^_ r ^-uv ( A (N-2\ M (N-1\\ t T%_ lN { , {N\ f'X- l 

v V ^v y _ 2V' r V Jv y ^ v Jv 77' 2 V Jv 7 ^ V Jv 

We recall that T 0,i (z) can be viewed as a probability measure on [0, 1] by setting 

JV 

r( 

N 



N J 2 \ \NJ^\NJJ 2 y \N J \ N 
(2N-2\ Y c 21 ( /2JV-3\ /2JV-2\\ T c 21 ( ( 2N — 1\ ,(2N-2 



3=1 



for all bounded continuous function / on [0,1]. Let < r <C 1 and / be a smooth function with 
support included in [0, 1 — r). In order to make the link between T°' l (z) and the process X N , let us 
introduce an extension of / by setting 

{f(-v+l/N) if v G [- (JV- l)/N,0) 

f(v + l/N) if wG [0,(JV-1)/JV] 

f(-v+{2N-l)/N) if v e [(JV- 1)/JV,2(JV- 1)/JV] 

/(w- (2JV-3)/JV) if »6 [2(JV-l)/JV,(3JV-3)/JV]. 



With these two functions we get the following representation. VZ G {1, . . . , JV}, 



T / °' i (z)=E^[/^(Xf)] 



Moreover, we have 

r/^z) = Em [r (xf )] + o (A) = [r ))] + o (i ; , 

where 



&•(«) = 

with v s G (1 — r, 1 — r/2), and where 



w if w G (—(1 — r), 1 — r) U (1 + r, 3 — r) 

u s elsewhere, 



/(-«) if ue[-i,o] 

/(«) if »e [o, l] 

/(-u + 2) if we [1,2] 

f(v-2) if we [2,3]. 



Let u G [0, 1) such that Z(JV)/JV — ► u. One can assume we [0, 1 — r) by changing r if necessary. As 
in the proof of Theorem |5.2[ we have the following lemma. 
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Lemma 6.14 Vip G Cg°(K). 



lim sup 

N^+oo velN 



L if 



£a r ,oo<P( v ) 



= 0, 



where 



In = 



u 



N — 1 — \Nr] 1 

-, U 

N ' N 

~N-l + [Nr] 2N-3 



N ' N 
and a rj00 is a ^-extended version o/ooo such that 



1 N - 1 - [JVr] 
iV' iV 

2JV - 1 37V - 3 - [iVr] 



u 



N ' 



N 



a r ,oc(v) 



aoo(-v) if «e(-(l-r),0] 

aoo(w) «/ u S [0,1 - r) 
aoo (-« + 2) »6(l + r,2] 

aco (w-2) »e[2,3-r), 



and i/ie martingale problem associated to £ ar x and starting from u is well-posed. 

Lemma 6.15 The law of the process (g r (X N ))^ starting from d(N) = (l(N) — 1)/N is tight on 
£>([0,+oc),R). 



Proof (of Lemma [6Tl5| ) Let T t N = cr(X? , s<t). According to Theorem 3 in [TBJ Chapter 3]. We 
have to show only the two following points. First, we have 



lim lim 



d(N) 



suv\g r (X t N )\>K) = 0, 

t>o 



since VJV, swp t>0 \g r (X^)\ < 3. Second, we have for each N, h € (0, 1), s € [0, h] and t > 0, 



E d{N) {(g r (X t N +s ) - g r (X t N ))W) < Kh 



In fact, we have 



E d{N) ((g r (Xy +s ) - 5r (Xf )) 2 |^f) < 2E d{N) ((M^t + s) - M^(i)) 2 |^ 



2E 



d(N) 



t+s 



£ N g r (X»)dw 



■N 



with 



M?(t) = g r {X?)-9r(XS')- C"9r{X?)d. 



which is a (T i r ) t >o-martingale. We also have 



since by Lemma 6.14 



sup sup \£ g r (v)\< +oo 



sup sup | < +oo. 



Moreover, £ N g r (0) = £ N g r (2(N - 1)/N) = 0. Then, we get 



E 



d(N) 



t+s 



£ N g r (X»)dw 



< Ch 2 



We recall that 



< Ml > t = f {£ N g r 2 ~ 2g r £ N g r ) (Xf)ds. 
Jo 
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Consequently, by the martingale property of (M^(t)) t >o, 

E d(N) ((M^(t + s) - M»(t)) 2 \F t N ) = E d(JV) ((M»(t + s) - MZ{t)) 2 \T») 

= E d{N) (M£(t + s) 2 -M»(tm N ) 
= E d{N) (< M» > t+s - < Mf d > t 



rN 

9. 

r t+s 



= E 
< Ch. 

In fact, in addition to the previous arguments, we also have 

sup sup \L N gl{v)\ < +00, 
N vei N u{v s } 



d (N) (J S ( £N 9r - Wor) (X£ )dw 



N 



sup N C N g*(0) = j$ < +00, and sup^ £ N g?(2(N - l)/2) = 2S# < +00. That concludes the proof 



: N 9 2 r (Q) - ^ 
Lemma T6. 151 □ 

Now, let us introduce some notations. Vj E N*, let 

r r w = inf (t > t^-V, x(t) e [-1, -(1 - r)) U (1 - r, 1 + r) U (3 - r, 3]) 
r« = inf (t > t?\ x(t) e (-(1 - r), 1 - r) U (1 + r, 3 - r)), 

with Tr,c = 0. Using the previous lemma, there exists (N') such that 



lim E d(NI) [f s (g r (X? ))] = E^ [/' (*(*))] . 

JV' — »-}-00 



Moreover, 



eQ" [r( X (z))] = £eq- [n^i^^a 

= ^E^[E^[r(,(z))l (T ,_ 1) ^ <T « 3)) 



3>1 



where A'Jt = ct(x(s), < s < t). With the following lemma we can identify each excursion between 
r r ( f- X) and 4 j) . 

Lemma 6.16 Vj G N*, the conditional law Q u ( • |.M u-i)) coincide up to the stopping time Tr 
with the conditional law P u ( • | ^Vt 0-1)), where P^ is £/ie unique solution of the martingale problem 
associated to £ ar x and starting from u. 



follows from Lemma 
Lemma 



6.14 



6.14 



Proof (of Lemma 6.16 ) This proof is a conditional version of Lemma 6.12 Moreover, this lemma 

and the fact that we are studying excursions between t^ 1 ^ and Tr . By 



in addition to g r (X?) = X* for rife 1] < z < t, 



Jj) 



limE 

N 



d(N) 



C^X^-Car^X?) d S 



0, 



and we also have 



En 



l(x^=o)dt 



L (Xf=2(AT-l)/Ar)( 



=0 




1 






A(l- 




=0 




1 






A(l- 


«')) 



-V(N-l)/N 

by symmetry of the process J*. As in the proof of Lemma 6.12 we get that V</5 € Co°(R) 



tAr, 1 . 
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is a martingale under the conditional law Q u ( • l-M^-i)). Finally, from the uniqueness of the mar- 
tingale problem associated to L aroo , Q u ( ■ \M (j-i)) coincide up to the stopping time with 
P u ( ■ |A / l T (j-i)) (see Theorem 6.2.2 in [27] )■ That concludes the proof of Lemma 



6.16 



□ 



From the previous lemma, Vj € N*, we have 



/*(*(2))l (Tr W-i)<, <T 0)) 



= E ff 



/ s (x(z))l (r , J - 1) < z<T ( 3)) 



and then 



hm E d(NI) [f(g r (X? ))] = E^ [f(x(z))] = E p - [f s (x(z))] , 

iv — >+oo 



where the limit does not depend to (N'). Consequently, 



hm T^\z) = E^[f s (x(z))] =T f (z,u), 



with 



8 

dz 



T f (z,u) = C art JT f (z,u) = C aoa T f (z,u). 



For the boundary conditions, first let ft, € (0, 1) such that < h <C 1, we have 

i(T / (z,/ l )-r / ( 2 ,-fc)) = i hm (E[^[/ s (Xf)] -E_ m [/ s (Xf)]) = 0, 
because of the symmetry of the process X N and f s 7 and therefore, 

2^T /M ) = C, 

Second, in the same way, let ft, € (0, 1) such that ft <C 1. Moreover, one can assume r < ft by changing 
r if necessary. Then, we have 

l(T f (z, 1 - ft) - T f (z, 1 + ft)) = i ^lim^ (E [w(1 - fe)1 [/ s (Xf )] - [/ s (Xf )]) = 0, 



and therefore, 



8 



2-T /( ,,l) = 0. 



As a result, using the density of the smooth functions with compact support in L 2 (0, 1) for ||.||l 2 (o,i) 
and the dominated convergence theorem we get the first point of Theorem |5.4| The second point 
is a consequence of the maximum principle and the density for the sup norm over [0,1] in {cp € 
C°([0, 1]), if(l) = 0} of the smooth functions with compact support included in [0, 1). ■ 
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